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SUMMARY 

^PROP5  is  the  third  version  of  a  computer  program  for  refining  the 
orbital  parameters  of  an  earth  satellite.  A  description  of  the  program  is 
given,  together  with  full  instructions  on  how  to  use  it. 
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INTRODUCTION 


PROP  is  the  acronym  for  a  computer  program  for  the  refinement  of  the 

orbital  £arameters  of  an  artificial  earth  satellite.  It  vls  conceived  in 

12  3 

July  1965.  as  an  eventual  successor  to  the  Pegasus  progrsuns  '  which  were 
used  at  R.A.E.  from  I962  to  The  new  program  has  been  written  in  the 

FORTRAN  language  (A.S.A.  standard  version  )  so  that  it  may  be  used  on  most 
large  modem  computers.  It  was  developed  on  the  Atlas  computer  of  London 
University  and  is  now  working  on  the  ICL  1907  computer. 

PROP  is  based  on  Merson's  analytical  model^  for  the  orbital  motion  of 
a  satellite.  The  parameters  of  this  motion  are  refined  by  use  of  observations 
of  the  satellite,  and  observations  of  many  tjT)es  may  be  used,  such  as  range, 
right  ascension  and  declination,  direction  cosines,  rauige  rate,  etc.  The 
present  paper  describes  some  of  the  structural  features  of  the  program  and 
gives  detailed  instructions  for  using  it. 

In  a  preliminary  version  of  the  program  -  the  unpublished  PROP!  -  it 
was  only  possible  to  determine  a  single  set  of  orbital  parameters  from  a  single 
set  of  observations  of  one  satellite;  further  determinations  had  to  be  set  up 
ab  initio.  With  a  later  version  -  PR0P2  -  it  became  possible  to  deal,  in 
automatic  succession,  with  sets  of  observations  of  a  number  of  satellites.  In 
the  current  version  -  PR0P5  -  a  large  numb^jr  of  fairly  minor  improvements  has 
been  introduced. 

PR0P3  consists  of  a  main  program  and  58  subprograms,  as  listed  in 
Appendix  A.  The  structural  relationship  of  the  program  units  is  illustrated 
at  Fig.1  «md  a  flow  chart  for  the  main  program  is  provided  at  Fig. 2.  T’le 
program  is  available  to  interested  parties  either  through  the  FORTRAN  source 
listings  or  as  a  source  deck  of  cards  punched  in  the  Atlas/lBM  Fortran  code, 
also  known  as  the  ’’bed"  code.  A  complementary  program,  known  as  PREP,  is  also 
available  ;  this  is  a  program  for  reporting  ephemeris  position,  based  on  the 
PROP  dynamical  model  and  Incorporating  most  of  the  PROP  subprograms. 

This  Manual  may  be  read  straight  through,  but  is  intended  mainly  as  a 
work  of  reference.  To  this  end  an  index  has  been  provided.  Also,  comprehensive 
cross-references  have  been  Included;  thus  a  reference  to  section  5  or  section 
6.1  would  be  indicated,  respectively,  by  "(S3)''  or  "(S6.1)". 
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PROGRAM  FUNCTION 


The  function  of  PROP5  may  be  suirunarlzed  as  follows:  given  estlsiates  of  the 
orbital  elements  of  an  .  *th  satellite  at  some  specified  epoch,  and  given 
observations  of  the  satell  te  over  a  period  of  some  days,  the  program  refines 
the  estimates  of  some  or  all  of  these  elements  by  an  iterative  least-squares 
differential-correction  procedure.  The  Initial  elements  are  provided  either 
directly  by  the  user  or  else  by  the  program  itself,  in  the  latter  case  by 
prediction  from  elements  at  some  other  epoch.  By  means  of  this  prediction 
facility,  orbital  elements  at  a  series  of  epochs  can  be  generated  during  a 
single  visit  to  the  computer;  a  set  of  observations  is  provided  for  each  epoch 
of  the  series,  but  estimated  initial  elements  are  required  for  the  first  epoch 
only. 

2.1  Dynamical  model 

During  the  orbit  refinement  PROP  takes  account  of  perturbations  due  to 
the  harmonic  coefficients  of  the  earth's  gravitational  field  and  to  atmospheric 
drag.  Formulae  for  the  short-periodic,  secular  and  long-periodic  effects  of 
the  earth's  zonal  harmonics  have  been  given  by  Merson^.  The  effect  of  the 
tesseral  harmonic  J,  „  on  the  mean  anomaly  of  the  satellite  is  taken  into 

d.tL  2 

account  by  the  method  described  by  Gooding  . 

The  main  effect  of  atmospheric  drag  on  a  satellite  is  to  reduce  the  semi- 
major  axis,  a,  and  hence  to  increase  the  mean  motion,  n.  This  effect  may 

a 

be  allowed  for,  in  PROP,  by  taking  the  rate  of  change  of  mean  motion,  n,  as 
a  parameter  to  be  determined.  The  program  calculates  the  small  effects  of 
drag  on  the  elements  e.  i,  0  and  (eccentricity,  inclination,  right 
ascension  of  the  node  and  argument  of  perigee)  in  terms  of  n. 

No  account  is  taken  of  perturbations  due  to  the  sun  and  moon.  If,  how¬ 
ever,  the  secular  effects  of  these  can  be  calculated  beforehand,  constant 
contributions  to  the  rates  of  change  of  the  orbital  elements  can  be  fed  into 
PROP.  In  the  next  version  of  PROP  it  is  hoped  to  incorporate  such  calcula¬ 
tions  internally.  It  is  remarked  that  constant  rates  of  change  are  only 
valid  for  periods  up  to  about  three  days,  but  that  for  close-earth  satellites 
l\ul-solar  perturbations  are  very  small. 

2.2  Observations  and  obseived  quantities 

An  observation  consists  of  data  obtained  by  a  sensor  (ground-station) 
at  a  particular  time.  It  consists  of  one  or  more  observed  quantities,  e.g. 
range  rate  or  a  pair  of  direction  cosines.  The  differential-correction 


(ubprogrun  of  PROP  can  handle  the  following  16  type*  of  obaervatlon.  but  only 
Types  1.  2.  J,  ^  and  7  are  accepted  by  the  PR0P5  Input  subroutine. 


Obse.-vei  quantities 


1 

p  (range) 

2 

n,  f  (right  ascension  i 

5 

^  m  (west-east  and  soi 

u 

p  (rsjige  rate) 
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f  (R.A./dec.  rates) 
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^  m  (direction  cosine 
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P  0  ?■ 
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p.  m 

9 

•  •  • 

p.  a  f, 

•  •  « 

10 

p.  i .  IT 

11 

p.  P 

12 

c.  0.  i 

15 

t,  m.  1,  in 

U 

p.  0.  6,  p,  a,  i 

15 

p .  i,  m,  p,  1,  m 

16 

X.  y.  z.  X,  y.  z. 

Type  16  does  not  relate  to  observations  in  the  normal  sense;  it  has  been  Included 
In  PROP  for  compatibility  with  the  epheraerl* -generation  program  PREP^. 


Observations  are  read  from  cards,  each  of  which  is  punched  in  one  of  a 
number  of  possible  formats  and  the  program  allocates  the  appropriate  type 
number.  Azimuth  and  elevation  are  converted  to  right  ascension  and  declination, 
and  so  come  under  type  2.  7.  12  or 

The  observation  formats  are  associated  with  the  sources  of  the  observa¬ 
tions.  l.e.  the  organizations  concerned  rather  than  with  the  observation  types. 
Descriptions  of  the  current  fonr.ats  in  which  observations  can  be  punched  are 
given  in  Appendix  B.  The  input  subroutine  can  easily  be  extended  to  cover 
further  formats. 


2.3  Observation  weights 

The  differentlal-coirection  procedure  uses  the  weighted  least-squares 
technique.  Every  observed  quantity  is  weighted  according  to  an  a  priori 
estimate  of  its  accuracy  (one-sigma  level)  -  the  weight  is  in  fact  taken  equal 
to  the  inverse  square  of  this  estimate. 
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For  ctrtain  groups  of  sonsors  standard  weights  are  used.  Observations 
trcm  sensors  in  these  groups  are  automatically  weighted  from  fixed  values  of 
the  accuracy  sstiaates;  these  values  are  set  by  the  main  program.  The  observa¬ 
tions  currently  treated  in  this  way  are  as  follows 

(i)  Minitrack  observations  {t,  m);  the  estimated  accuracy  of  a 
direction  cosine  is  0.00029>  equivalent  to  an  angular  accuracy  of  1  min  arc. 

(ii)  NASA  range  and  range-rate  obaervationa  (p,  the  estimated 
accuracies  are  I50  m  and  0.^  m/s  respectively. 

(iii)  Moonwatch  visual  obaervationa,  as  received  on  punched  cards;  the 
estimated  accuracy  is  0*.1,  interpreted  aa  O.OOl?  radian. 

For  all  other  observations  the  estimate  of  accuracy  must  be  punched  for 
each  observation  as  part  of  the  format. 

When  more  than  one  observed  quantity  is  included  in  an  observation,  the 
quantities  are  regarded  as  Independent,  The  estimated  accuracy  of  both  ^ 

direction  cosines  of  a  pair  is  taken  as  the  same.  With  Q  and  t,  however,  the 
estimated  accuracy  refers  to  6 ;  the  estimated  accuracy  of  o  is  given  by 
0(0)  =  0(6)  /  cos  i.  It  is  because  this  weighting  of  right  ascension  and 
declination  is  equivalent  to  an  analogous  weighting  of  azimuth  and  elevation 
that  azimuth/elevation  observations  are  converted  to  right  ascension/declina¬ 
tion  and  not  given  their  own  type  number. 

Time  is  not  treated  as  an  observed  quantity  in  PROP,  but  for  observations 
which  have  a  priori  estimates  of  both  timing  and  angular  errors  (0^,  o^),  the 
estimated  accuracy  in  &  is  sissle,iied  a  compounded  value  given  by: 

0(6)  =  [0^  +  J  &  0^^  , 

a 

where  0  is  an  approximation  to  the  mean  tracking  rate,  given  by 

=  n  /  (1  +  -  R/a)  , 

# 

where  n,  e  and  a  are  mean  motion,  eccentricity  and  reml-major  axis  of  the  orbit 
and  R  is  the  mean  equatorial  radius  of  the  earth. 

2.^»  Sensors 

The  geographical  location  of  sensors  must  be  provided  on  punched  cards, 
one  card  for  each  sensor.  Instructions  for  punching  sensor  cards  are  given 
in  section  0,6, 


J 
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3  TIME  AWD  COORDIMATE  8T8TEK3 

Th«  bade  time  lystem  uaed  by  PROP  is  UTl  and  epoch*  should  normally 
be  understood  to  be  UTI  epochs.  The  times  of  many  observations  are  given  In 
the  UTI  system,  but  some  observations  (e.g.  Mlnltrack  observations)  refer 
to  UTC,  the  system  defined  by  WWV  transmissions,  and  others  (e.g.  Baker- 
Miiwn  camera  observations)  to  the  atomic-time  system  A1 .  The  input  routines 
include  the  facility  to  correct  WWV  or  A1  times  to  ITTI  .* 

Calendar  dates  are  reckoned  in  Modified  Julian  Days  (MJD)’.  Taese  are 
related  to  (ordinary)  Jhllan  Days  by  the  formula  MJD  ^  jp  _  2400000,5. 

Thus  during  the  present  century  each  midnight  is  represented  by  a  fire-digit 
integer;  e.g.  i960  JAN  1.0  is  represented  by  36954. 

PROP  uses  a  system  of  rectangular  coordinates  defined  as  follows.  The 
origin  0  is  at  the  centre  of  the  earth  and  the  axlt'  Oz  points  towards  the 
north  pole.  Ox  lies  in  the  plane  of  the  equator  of  date,  but  instead  of 
pointing  towards  the  true  equinox  of  date  it  points  towards  a  projection  of 
the  mean  equinox  of  the  epoch  1930.0.  More  precisely,  the  position  of  Ox  is 
obtained  by  rotation  from  the  true  equinox,  by  the  amount  of  the  precession 
and  nutation  in  R.A.  since  1930.0,  but  in  the  opposite  sense.  Oy  completes 
the  right-handed  system  Oxyz. 

The  above  system  is  the  one  used  by  the  Smithsonian  Institution 

0 

Astrophyrlcal  Observatory  in  the  D.O.I.  program  .  It  is  remarked  that  1930.0 
above  actually  means  22h  09m  42s  on  1949  EEC  51*  i.e.  55281.9254  MJD. 

For  further  Information  on  systems  of  time  and  coordinates,  in  relation 

9 

to  satellite  orbits,  the  comprehensive  treatment  by  Veis  should  be 
consulted . 

4  ORBITAL  ELEMENTS  AND  ORBITAL  PARAMETERS 

4.1  Orbital  elements 

In  the  orbital  model  used  in  PROP  there  are  five  basic  elements: 
eccentricity,  e;  inclination,  i;  right  ascension  of  the  node,  0; 
argument  of  perigee,  and  mean  anomaly,  M.  As  In  the  D.O.I.  program 

these  are  mean  elements  in  the  sense  that  they  are  free  of  first-order  short- 
periodic  perturbations.  Variation  of  the  basic  elements  is  represented  by 
polynomials  in  time,  with  additional  perturbation  terms;  thus 

*In  orbits  determined  from  Minitrack  observations  this  facility  has  not 
so  far  been  used;  i.e.  orbital  parameters  refer  to  the  UIC  system  in  which, 
it  should  be  noted,  there  are  sometimes  IOO  msec  "jumps". 
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e  ■  e  4  de  , 

1  -  I  ♦  dl  , 

n  =  ff  4  dO 

CJ  =  .'  4  <>*' 

and  M  c  M  4  dM  . 

vhare  e  etc.  are  polynomials  in  time.  i.e. 

e  »  V  t*^  .  etc. 

Here  e  .  i  .  0  .  .  M  are  values  of  the  mean  elements  at  epoch  (when 

O  O  0  o  o 

t  o)  and  e.  i.  d.  (j,  m  are  values  at  time  t.  i.e.  the  perturbations  de 
etc.  are  zero  at  epoch,  (it  is  remarked  that  e^.  i^  etc.  are  not  in  general 
the  values  of  the  osculating  eccentricity,  inclination  etc.  at  epoch,  since 
the  first-order  short-periodic  perturbations  in  the  osculating  elements  are 
not  defined  to  be  zero  at  epoch  -  see  also  Appendix  C.)  The  de  etc.  do  not 
appear  explicitly  to  the  PROP  user,  who  need  not  normally  be  concerned  with 
them.  In  PROPJ  they  represent  long-periodic  perturbations  due  to  the  zonal 
harmonics  of  the  earth's  gravitational  field  and  atmospheric  drag,  but  they  may 
be  extended  to  cover  the  tesseral  harmonics  and  luni-solsr  effects,  in  a  later 
version  of  the  program. 

The  degree  of  each  polynomial  is  optional,  between  the  values  0  and  ^ 
(inclusive),  except  that  the  lower  limit  of  the  degree  of  M  is  1,  since 
H|  -  the  mean  motion  at  epoch  -  must  always  appear.  The  polynomials  provide 
for  the  secular  variation  of  the  elements  due  to  the  earth's  zonal  harmonics 

and  atmospheric  drag.  The  set  of  coefficients  e^,  e^,  ....  i^,  . etc. 

in  the  dynamical  model,  between  7  and  30  of  them  altogether,  will  be  referred 
to  as  "the  orbital  elements"  of  the  satellite. 

The  mean  value  of  the  mean  motion,  n.  is  a  derived  parameter  given  by 
the  rate  of  change  of  M;  thus 

'  Z  • 


n 


where  -  (j  +  1 )  Mj^, . 

The  leBl-BeJor  axle  e  li  defined  by 


•  ■  -  i  Jg  H*  (M/n®)‘j  (1  -  U  iln^IXl  - 


where  m  ie  the  eerth*a  greyltetional  conitent,  le  Its  second  tonal 

hareonic  coefficient  and  R  is  its  mean  equatorial  radius. 


4,2  Orbital  paraaeters 

The  complete  set  of  orbital  elesMnts,  as  Just  defined,  are  sosMtiaes 
referred  to  as  the  orbital  paraaeters  of  a  satellite.  With  PROP,  however,  it 
is  convenient  to  reserve  the  word  "parasMters*  for  a  subset  of  the  set  of 
eleaents,  naaely.  those  elesMnts  which  are  refined  by  the  differential- 
correction  procedure.  At  least  one  eleaent  aust  be  a  paraaeter  -  since  othar- 
wlse  there  would  be  nothing  for  PROP  to  do  -  and  the  najclaua  number  of 
parameters  allowed  has  arbitrarily  been  set  at  20.  thoufh  it  would  be  unueual 
to  use  more  than  1o.  It  la  possible  for  all  eleaents  to  be  paraasters.  Any 
that  are  not  are  regarded  as  being  *held  fixed*,  whereas  the  paiaaetere 
aay  be  thought  of  as  "floating*. 


say 


The  situation  aay  be  clarified  by  consideration  of  the  e 
K-1 

e  •  V  e^  where  1  *•  K  !■  6,  in  which  caee  there  are 


polynoalal. 

%  eleaenta 


associated  with  e.  Then  the  first  n.  for  any  a  such  that  0  ’  a  ^  L 
aay  be  chosen  by  the  PROP  user  to  be  paraaeters  and  the  reaalnder.  if  any,  to 
be  "fixed  eleaents".  It  is  not  possible  to  have  a  paraaeter  with  suffix 
higher  than  that  for  any  of  the  fixed  eleaents.  e.g.  to  have  a|  as  a  para¬ 
aeter  while  bolding  e  fixed. 

o 

At  the  beginning  of  each  PROP^  run  a  coaplete  set  of  values  of  the  orbital 
eleaents  is  sMde  available  to  the  coaputer.  and  the  parameters  are  refined  at 
the  end  of  each  iteration.  The  "fixed*  elements  are  in  general  held  fixed 
throughout  the  run.  as  the  name  suggests.  In  the  case  of  e^,  and/or 

however,  contributions  to  these  elements  are  coeiputed  afresh  at  the 
beginning  of  each  Iteration,  froa  foraulae  which  represent  the  theoretical 
effects  of  certain  perturbations,  so  these  elements  will  vary  slightly  froa 
one  iteration  to  the  next  (S4.3). 
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Althoui^  there  li  coniiderable  flexibility  In  the  possible  choice  of 
elffMnts  end  pareiteters^  there  ere  two  main  ways  In  which  PROP?  will  be  used 
and  those  are  now  eusmtarlzed. 

Standard  B)odel  -  observations  are  available  over  a  period  of  a  few  days, 
ueually  not  B»re  than  ^  or  5:  e.  T.  n  and  »-»  are  linear,  while  M  Is  quad¬ 
ratic  or  cubic;  e  ,  1  ,  n  ,  ,  m  ,  M,.  Mr,  and  M,  (if  R  is  cubic)  are 

the  parameters  with  M2  (and  M^)  providing  an  empirical  representation  of 
acceleration  due  to  air  drag;  e^.  1,^.  ere  "fixed"  elements; 

there  are  11  (or  12)  elements,  7  (or  ^)  parameters. 

Alternative  model  -  observations  over  a  much  longer  period,  perhaps  more 
than  a  month  ere  to  be  combined:  the  {>olynomlals  are  of  high  enough  order  to 
fit  the  observations,  and  al 1  the  elements  are  parameters, 

^ . 5  Exclusive  and  inclusive  elements 

The  elements  e,^  1,^  and  .j  will  almost  always  be  present  in  the 

dynaaical  model  chosen  for  the  orbit.  They  represent  the  secular  rates  of 

change  of  e,  1.  fl  and  •,  and  they  receive  special  treatment  by  PROP, 

Consider  the  e  polynomial,  say,  Whereas  e^  and  such  of  e^,  e^.  e^  and 
as  are  in  the  dynaaical  model  are  either  treated  as  parameters  or  else 
held  absolutely  fixed  at  their  initial  values,  for  e.  the  situation  Is  sore 
cosplicated.  This  element  as  used  to  evaluat#*  the  polynomial  as  a  function 
of  time,  is  forced  by  the  addition  of  two  components.  One  of  these  components 

has  the  saae  status  as  e  .  e.  etc.  in  that  an  arbitrary  value  may  be  read 

o  < 

Into  the  computer  initially  and  then  treated  as  a  parameter  or  held  fixed. 

This  component  Is  called  the  "exclusive"  e.^  ;  it  will  normally  be  much  the 

smaller  of  the  two  components  and  Indeed  will  frequently  -  like  Cp  etc.  - 

be  zero.  The  other  component  is  compute!  as  a  function  of  the  elements  e^, 

1  .  w  M,  and  (assuming  M,  present  in  the  model),  at  the  beginning  of 
o  O  1  t  C 

every  iteration  of  the  PROP  differential-correction  procedure.  This  component 
has  not  been  given  its  own  name  but  its  combination  with  the  exclusive  e^ 
is  called  the  "inclusive"  e.^  .  For  consistency  of  definition  it  has  been 
arranged  that  the  cor.blnatlon  of  a  '-omputed  coteponent  of  e^  with  the 
exclusive  component  is  always  carried  out  (assuming  e.^  is  present  in  the 
model),  even  when  e^  is  being  treated  as  a  parameter. 

The  above  explanation  for  e.^  applies  also  to  i^  0.^  and  .  The 
computed  coftponents  of  these  elements  allow,  in  PROP?,  for  secular  perturba¬ 
tions  due  to  the  earth'i  zonal  harmonics  and  atmospheric  drag.  The  exclusive 
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elements  represent  other  secular  perturbations;  they  may  be  treated  as  para¬ 
meters,  or  set  to  zero,  or  given  fixed  values  obtained  by  the  user 
independently  of  PROP?  -  for  example  to  represent  luni-solar  secular  effects. 

The  values  of  e,^  etc.  on  punched  cards,  both  as  input  at  the 
beginning  of  a  run  and  as  the  final  output  provided  by  the  computer,  are 
exclusive  values.  The  values  given  by  the  lineprinter  output  (S9)  are 
inclusive  values . 

4.4  Further  remarks  on  exclusive  and  inclusive  elements 

Two  further  points  must  be  mentioned  in  order  to  provide  a  complete 
explanation  of  the  system  of  elements  used  in  PROP.  They  are  rather  technical 
and  the  rest  of  this  section  may  be  omitted,  without  loss,  by  the  general 
reader. 


First,  it  was  decided  at  a  late  stage  in  the  development  of  PROPJ  to 
extend  the  distinction  between  exclusive  and  inclusive  elements  to  e^  etc. 
with  j  >  1.  the  object  being  to  inrorove  the  representation  of  drag.  Prior 
to  this  modification  the  situation  was  simply  that  the  drag  sub-components  of 
the  computed  components  of  e.^  etc.  contained  the  acceleration  element  Mg 
as  a  factor;  i.e.  they  were  proportional  to  the  value  of  n  at  epoch.  But  a 
more  natural  procedure  would  have  been  to  arrange  for  drag  components  of  e 
etc.  to  be  proportional  to  n  all  the  time,  i.e.  to  introduce  computed 
components  of  eg  etc.,  proportional  to  ng.  and  so  on  for  as  many  n*s 
as  are  in  the  model.  Thus  the  modification  that  was  made  was  to  introduce 

exclusive-to- inclusive  computation  for  e  etc.,  for  as  many  j  as  there 

J 

are  both  e^  and  n^  in  the  chosen  model.  The  zonal-harmonics  sub-components 
still  apply  to  j  =  1  only,  and  no  provision  has  been  introduced  for  present¬ 
ing  the  computer  with  fixed  exclusive  values  of  e^  etc.,  with  j  ^  1. 
other  than  zero.  So  long  as  none  of  e.  i.  ^  and  i>)  is  represented  by  a 
quadratic  or  higher-degree  polynomial,  e.g.  so  long  as  the  standard  model  is 
used,  PROP  will  be  unaffected  by  this  modification. 


Second,  for  reasons  explained  by  Ctooding  ,  it  would  be  unsatisfactory  to 
permit  the  existence  of  exclusive  elements  0^  and  ( j  ?  1 )  in  the  orbital 
model  unless  certain  suijustments  made  to  the  elements  .  If  the  Mj 
presented  to  the  computer  are  regarded  as  "exclusive",  in  a  further  extension 
of  the  meaning  given  in  section  4.3,  then  the  adjustment  to  each  Mj  consists 
of  the  subtraction  of  Wj  and  Oj  cos  i^.  the  adjusted  Mj  being  regarded  as 
"inclusive".  Since  this  adjustment  must  apply  to  the  mean  anomaly  but  not 


n 


299 


to  the  mean  motion,  quantities  nj.  equal  to  (j  +  1)  are  computed  and 

stored  while  the  are  still  In  exclusive  form.  The  and  are  normally 

very  small  amd  hence  there  is  usually  not  much  difference  between  the  exclusive 
and  Inclusive  values  of  the  Mj  • 

It  Is  again  stressed  that  the  elements  printed  by  PROP  are  inclusive,  but 
that  elements  on  cards  are  always  exclusive. 

5  URITS  AND  CONSTANTS 

Distances  are  measured  in  kilometres.  Angles  are  expressed  in  degrees  for 
input  and  output,  but  are  held  as  radians  inside  the  computer  (though  this  is 
normally  of  no  consequence  to  the  user).  Time  is  in  days  for  the  input  and  out¬ 
put  of  orbital  elements  (e.g.  is  in  degrees/day  and  M2  in  degrees/ 

day/day),  but  otherwise  in  seconds  (e.g.  in  the  storage  of  orbital  elements 
and  the  output  of  range-rate  residuals). 

Care  has  been  taken  over  the  specification  of  constants  used  in  the  prograsi. 
The  need  to  chMge  the  values  of  the  geophysical  and  certain  program  constants 
has  been  anticipated  and  the  current  value  of  each  constant  is  assigned  once  only, 
either  in  the  main  program  or  else  in  the  block  data  subprogram:  communication 
between  program  units*  is  by  the  FOPTRAW  variables  to  which  these  values  have 
been  assigned.  To  assign  a  new  value  to  a  constant  it  is  normally  necessary  to 
recompile  the  main  program,  or  block  data  subprogram. 

5.1  Geophysical  constants 

current  values  of  the  geophysical  constants  are  assigned  in  the  block  data 
subprogram  an'’  listed  In  every  PROP  output.  Values  of  M.  J2  R  (S^.l)  are 
given  -  in  terms  of  the  FORTRAN  variables  used  -  as  follows  :- 

EMU  •  5986'2.^  (kinVs«c^)  EJ(2)  '  1  x  1  and  ERAD  •  6576.165  (ka). 

Values  of  J,  J. .  J.  r  have  been  assigned  from  the  results  of  Klng-Hele 

1)11*  ^  ^ 

et  al.  to  Ej(5)  etc.  The  maxlm.’T".  degree  of  the  tonal  harmonics  allowed  for. 

without  recoeipillng.  is  given  by  L  »  16.  (To  use  a  value  L  >  16  it  is 

BAX 

necessary  to  change  the  dimensions  of  the  arrays  EJ.  In  common  block  /ORBZT/. 
and  ABCD,  in  common  block  /PRECON/.  EJ  must  be  given  dimension  L  „  instead 
of  16  and  ABCD  swst  be  given  a  first  dimension  of  -  2  instead  of  14.) 

The  density  scale  height  of  the  earth's  atmosphere  is  currently  set  as 
DOfSCH  -25.0  (km).  The  values  of  the  tesseral  harmonic  J2  2 

•i  'procjwa  unli*  I*  eiilwr  Uw  mnin  pntcrno  or  or*  or  irw  ihrtt  cIasmi  of  mmtrocraai  'amroulirm*, 
■nmcilon*  WKl  'block 
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longitude  phase  used  in  the  mean -anomaly  correction  ,  are  given  by 

^  I  ^ 

EJ22  =  0.0000018  and  ELA1C2  =  -0.51  (radian)  reipectlvely. 

5.2  Program  design  constants 

The  Initial  value  of  c  (s6.6)  Is  normally  (S5.5)  ■et.  In  the  main 
program,  by:-  EPSLON  =  200.0. 

The  maximum  number  of  observations  (MAXIM)  and  the  maximum  number  of 
sensors  (MAXSTA)  which  are  allowed  are  set  to  8o  and  5"^  respectively.  If  it 
is  ae.'-irei  to  exceed  these  limits  -  or  to  reduce  them  to  save  storage  -  the 
main  program  must  be  recompiled;  the  statements  MAXIM  =  S'*  and  MAXSTA  =  30 
must  be  altered  and  so  must  the  Items  DATA(21  8o)  and  STASHN(0.  50)  in  the 
DIMENSION  list. 

Three  accept/reject  constants  are  set  by:  REJLEV(1 )  5.0. 

REJLEV(2)  -  U.-).  REJLEV(5)  =  l’'.^  They  have  the  function  of  dividing  the 
observations,  during  each  iteration  of  the  differential-correction  procedure, 
into  four  groups  (S6.7). 

5 . 5  Observation  weights 

The  weights  of  NASA  Minitrack,  range  and  range-rate  observations  and  of 
certain  optical  ^b.':o^\-at  i  ons  arc  given  standard  values  (S2.5).  The 
assigning  statements  are:-  SIGKA(1 )  -  0.00029.  SIGMA(2)  ^  0.I5. 

SIGMA(5)  -  '\0005  and  SIGMA('*)  ■  ^.I'l?. 

5 Input/output  constants 

The  card-reader,  lineprinter  and  card-punch  unit  numbers  are  set  In  the 
block  data  subprogram.  A  further  setting  gives  the  number  of  lines  on  a  line- 
printer  sheet,  l.e.  the  maximum  possible  number  of  lines  of  printing  which 
can  occur  assuming  no  top  or  bottom  margin.  For  running  on  Atlas  these  four 
Input/output  constants  are  set  by:  IF  ,  IP  X  IPUNCH  15  and 
ILINES  -51.  On  the  197  they  are  set  by:  IF  ■  1,  IP  *  2.  IPUNCH  -  5  and 
ILINES  -  Every  input  or  output  instruction  in  every  progra.m  unit  of 

PF0P5  refers  to  its  peripheral  device  by  the  appropriate  one  of  these 
FOFTRAN  variables. 

5 . 5  Constants  that  may  be  changed  without  recompiling 

New  values  of  four  of  the  constants  mentioned  earlier  can  be  set  with¬ 
out  recompiling  (S8.2).  The  four  are:  L  (must  be  >  16)  EJ22  DENSCH  and 
EPSLON  (Initial  value). 
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6  PyiR/gj,  PROGRAM  DESIQW 

6.1  Meanings  of  "run*  and  "Job* 

A  "run"  consists  of  the  complete  action  of  the  program  In  refining  orbital 
parameters  for  a  particular  satellite  at  a  particular  epoch;  the  epoch  is 
specified  by  an  integral  number  of  days  -  held  by  the  FORTRAN  variable  MJDOCH  - 
so  that  epochs  are  always  at  nh  (midnight).  A  run  starts  with  the  reading  of  a 
control  card  (S8.2)  and,  if  successful,  ends  with  the  output  of  the  final  values 
of  the  orbital  elements  and  of  the  covariance  matrix  of  the  parauneters.  How¬ 
ever,  a  run  may  end  in  failure  due  to  one  of  a  number  of  specified  causes;  for 
example  the  refinement  process  may  not  converge,  or  too  many  observations  may 
be  rejected  at  some  stage.  The  cause  of  failure  is  printed,  whenever  this 
occurs.  Besides  success  and  one  of  the  standard  failures  there  is  a  further 
possibility.  A  run  may  break  down  without  coming  to  an  "end"  at  all;  the 
most  likely  reason  is  that  the  data  deck  for  the  run  has  been  prepared 
incorrectly. 

A  "job"  consists  of  a  sequence  of  runs  carried  out  by  PR0P3  one  after  the 
other,  without  manual  intervention. 

6.2  Initial  elements  and  prediction 

For  ea''h  run  of  a  Job  PR0P5  requires  a  set  of  initial  orbital  elements, 
including  the  estimated  values  of  the  parameters  which  are  to  be  refined.  These 
elements  may  be  provided  in  one  of  two  ways:  either  by  the  user  in  the  data 
deck  for  the  run.  or  else  by  the  program  itself,  making  use  of  the  results  of 
the  previous  run.  In  the  latter  case  it  must  be  true  that  (a)  thie  run  is  not 
the  first  of  the  Job  and  (b)  the  previous  run  was  for  the  same  satellite. 

The  situation  Just  referred  to  involves  the  idea  of  "prediction"^.  Given 
values  of  the  n  ‘a’,  elements  of  a  satellite  at  one  epoch.  PROP  can  be  staked  to 
predict  their  values  at  another  epoch,  using  the  normal  formulae  for  orbital 
perturbations.  Predictions  can  be  either  forwards  to  a  later  epoch  or  back¬ 
wards  to  an  earlier  one. 

The  validity  of  the  prediction  formulae  is  limited  by  the  difficulty  in 
accounting  for  the  long-term  effects  of  drag.  Prediction  over  a  period  of  a 
week  is  usually  satisfactory  even  for  low-perigee  satillltes.  since  it  is  only 
initial  elements  which  are  predicted.  Prediction  over  a  number  of  weeks  should 
normally  be  avoided,  since  the  error  in  the  (along-track)  element  in 

particular,  may  be  so  large  sis  to  cause  the  refinement  process  to  diverge. 


J 


The  us\ial  application  of  prediction  is  when  initial  elements  for  a  run 
are  predicted  automatically  by  PROP  from  the  final  elements  of  the  previous 
nin.  A  second  application  arises  when  the  user  provides  PROP  vit^  an  initial 
set  of  elements,  but  wishes  the  analysis  to  be  based  on  a  different  epoch 
from  the  one  associated  with  these  elements.  This  happens  when  a  series  of 
runs  is  spread  over  more  than  one  Job  and  the  results  from  the  last  run  of  a 
Job  are  used  in  preparing  the  data  for  the  first  run  of  a  new  Job. 

The  various  possibilities  are  covered  by  the  value  of  the  FORTRAN 
variable  NEWSAT,  read  from  the  control  card  at  the  beginning  of  each  run  (S0.2). 

6.5  End  of  Job 

The  most  satisfactory  end  of  a  Job  occurs  when  all  its  runs  have  been 
dealt  with  successfully,  i.e..  a  final  (converged)  set  of  orbital  parameters 
has  been  printed,  together  with  standard  deviations,  at  the  end  of  each  run. 

After  the  successful  completion  of  any  run.  the  next,  if  there  is  one. 
is  automatically  started.  Otherwise,  the  action  of  the  computer  depends  on 
whether  the  run  came  to  a  standard  failure,  or  "broke  down" (S6.1). 

If  a  run  comes  to  a  standard  failure  then  the  control  card  for  the  next 
run  is  read.  If  this  control  card  specifies  the  prediction  of  elements  from 
the  run  that  failed,  then  the  new  run  does  not  proceed  beyond  this  point;  any 
further  runs  are  abandoned  and  the  Job  ends.  If.  however  the  new  run  supplies 
its  own  elements,  then  the  failure  of  the  preceding  run  is  irrelevant  and  the 
Job  proceeds.  Thus  the  success  of  the  last  run  of  a  Job  does  not  imply  the 
success  of  all  preceding  runs.  The  occurrence  of  a  standard  failure  is  marked 
by  a  printed  statement  about  the  failure  {S9.^)  and  this  always  starts  with 
the  word  "nir.C ONTItiUK". 

If  a  run  breaks  down  the  Job  is  abandoned  at  once,  no  further  run  being 
started.  Here  "breaking  down"  covers  two  possibilities.  The  run  may  fail  in 
a  way  which  is  regarded  as  non-standard,  but  yet  which  hu  been  provided  for 
by  the  program,  since  a  suitable  diagnostic  is  printed;  e.g.  "ORBITAL 
EUMEHTS  ARE  FOR  WRONG  EPOCH".  Alternatively,  a  completely  unanticipated 
failure  may  be  reached;  it  is  hoped  this  will  be  a  rare  occurrence. 

6.4  "Mode"  of  refinement 

If  the  initial  estimates  of  the  orbital  elements  are  good,  and  if  the 
coverage  of  the  observations  (i.e.  their  distribution  around  the  orbit  and  in 
time)  is  adequate,  it  should  be  all  right  to  start  refining  all  the 
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paruaeters  together  straight  away.  When  the  eatlmatei  are  poor,  hovcvei.  it 
Is  better  to  start  by  refining  one  or  two  of  the  paraaeters  only,  temporarily 
holding  the  others  as  fixed  elements;  the  latter  can  be  introduced  after  one 
or  more  iterations  with  the  restricted  model.  Two  parameters  which  should 
be  refined  from  the  beginning  are  and  M^.  in  order  to  fit  the  mean- 
anomaly  variation,  i.e.  to  get  the  "timing”  right. 

The  way  in  which  the  parameter  model  is  restricted  in  PROP  is  by  means 
of  what  is  called  the  "mode"  facility.  There  are  three  possible  modes  of 
operation,  amd  the  initial  mode  is  Indicated  by  the  value  of  the  FORTRAN 
variable  MODE,  read  from  the  control  card  for  the  run  (S8.2). 

MODE  ■-  ^  gives  the  mode  in  which  all  parameters  are  refined  immediately. 
MODE  "  1  gives  a  mode  in  which  at  most  the  four  parameters  e^,  M^'  and 
M2  are  refined  -  "at  most",  because  one  or  more  of  these  elements  may  not  be 
parameters  to  be  refined  at  all.  MODE  2  gives  a  mode  in  which  at  most 
M^  and  are  refined.  For  MODE  >  2  the  mode  is  the  same  as  for  MODE  ^  2. 
but  the  process  will  stay  longer  in  this  mode. 

At  the  end  of  each  Iteration  if  no  obsexvations  have  been  rejected 
(S^,7)  and  if  MODE  >  the  value  of  MODE  is  reduced  by  unity:  otherwise, 
the  next  iteration  is  carried  out  in  the  same  mode.  Since  a  xnui  is  only 
permitted  to  complete  its  convergence  in  the  zero  mode,  all  successful  runs 
will  eventually  reach  this  mode,  with  all  restrictions  on  the  parameter  model 
xremo'/ed . 

6,5  Singularities 

Steps  have  beer,  taken  to  minimise  the  effects  of  the  well-known 

2 

singularities  in  the  orbital  elements  -  tho.se  at  sin  i  =  ,8  (the  critical 
inclinations),  at  e  0  (circular  orbits),  and  at  i  =  0*  or  l8o* 

(equatorial  orbits).  The  formulae  of  Merson^  for  the  position  and  velocity 
of  a  satellite  at  arbitrary  time  have  been  so  expressed  that  none  of  these 
singularities  causes  any  trouble  in  the  par^  of  PROP  which  computes  position 
and  velocity.  It  should  be  remarked,  though,  that  the  orbital  model  may 
represent  the  pex^urbations  in  the  satellite  motion  less  accurately  near  the 

12  n  H 

singularities.  Thus  Jook  has  shown  that  near  the  circular  singularity 
perturbations  are  more  accurately  expressed  in  terms  of  elements  e  cos  w 
and  e  sin  than  in  terms  of  the  normal  e  and  w. 

There  are  real  difficulties  associated  with  the  differential-correction 
part  of  PROP,  however,  at  both  the  "cl  ar"  emd  the  "equatorial"  singularities. 


These  IlfflcuUles  arise  from  the  fact  that  ^  and  M  became  indeterminate  at 
the  circular  si n<?u lari ty.  and  n  and  at  the  equatorial  one.  The  effecta 
are  t«_-fold:  first  the  convergence  of  the  process  may  be  adveraely  affected  - 
It  may  break  down  altogether  -  due  to  the  use  of  inappropriate  elementa  near 
a  singularity:  second  the  accuracy  of  the  computer's  basic  matrix-lnveralon 
operation  may  be  affected  by  the  fact  that  partial-derivative  matrices  become 
Ill-conditioned  due  to  the  presence  of  almost-equal  columns.  Either  or  both 
of  these  effects  may  be  minimised  If  the  PROP  user  gives  an  appropriate  value 
to  the  control  parameter  (FOPTRAH  variable)  JELTYP  (S“,2), 

Satisfactory  convergence  is  restored  by  a  temporarj'  transformation  of 

orbital  ;  orometers  at  .^ust  one  part  of  the  program..  Consider,  for  example,  a 

nearlv  equatorial  orbit  for  which  both  1  and  T  are  (floating)  parameters 
'  o  0 

for  wMich  corrections  .'1  and  have  u.-.t  been  computed.  Then  If 

0  0 

JELTYP  has  be^n  appropriate ly  set  th**  computer  does  not  correct  1^  and  0^ 

directly,  but  forr.s  corrections  to  the  temporary  param.eterr. 

sin  1  cos  and  sin  1  sin  H  In.stead.  using 
0  0  0  o 

in  1  co.c  )  C. i  -'OS  i  ''OS  -  .'il  sin  1  sin  ^ 

0  O  O  O  00  0  0 


an  i 


'  (sir. 


;  in 


.^'1  cos  1  sin  r  sin  1  cos  P 

o  o  0  0  0  0 


It  then  computes  n-'w  values  of  1  and  C  from  the  corrected  values  of  the 

(JO 

temporary'  para.meters.  Simllorlc.'  for  a  nearly  circular  orbit  the  parametera 


e  cos  e  sin  •  and 
0  o  o  o  0 


M  terporari  Iv  replace  e  and  M  .  For  an 

o  o  o  o 


orbit  which  is  both  nearly  equatorial  and  nearly  circular  the  parameters 


an' 


or»*  te'"’.  eras'  r'’:  a'^i  b'; 


-o  'c  '’o  *  V  ^o  ^  "’o  '  “^0'  ''‘’*>** 

k  is  a  fixed  nu-mber  equal  to  cos  1  and  it  Is  assumed  that  e  etc  are 

o  o 

Indeed  all  parameters  under  refinement. 


The  other  aspect  of  the  singularities  ill-conditioning  Is  also 

controlled  by  a  transformation  of  parameters.  Suppose  that  w  and  M  are  both 

o  o 

parameters  and  that  we  are  near  the  -ir  ulnr  singularity.  Then  partial 
derivatives  with  respect  to  these  parameters  are  almost  equal.  However  if 
these  parameter.-  are  replaced  by  ^  and  M  •  the  derivatives  with  respect 
to  became  v<*r/  small  and  the  matrix  of  partial  derivatives  is  no  longer 
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Ill-conditioned.  The  computer  now  obtains  correctlone  tow  and  m  w  from 

0  0  0 

which  the  correction  to  Is  of  course  Inmedlate.  Similarly,  near  the 

equatorial  singularity  derivatives  with  respect  to  and  are  replaced  by 

derivatives  with  respect  to  0  and  w  -«  k  0  ,  where  k  is  the  same  as  In 

o  o  o 

the  last  paragraph.  (Remember  -  a  partial  derivative  with  respect  to 
keeping  fixed,  Is  not  the  same  as  a  partial  derivative  with  respect  to 

n^,  keeping  ■*  k  fixed'  )  Finally,  when  near  to  both  singularities 
the  trans  format  ion  is  from  n,w.M  to  n.u  4kn,M 
It  may  be  objected  that  it  does  not  help  to  replace  an  ill-conditioned 
matrix  by  a  well-conditioned  matrix  containing  a  nearly  zero  column,  since 
matrix  inversion  will  then  lead  to  floating-point  overflow.  This  only  occurs 
when  e^  almost  exactly  equal  to  zero,  however,  and  in  this  remote 

contingency  PROP  will  fall. 

Correction  for  both  the  above  aspeci.s  of  the  singularities  is  by  setting 
JELTYP  (s8.2).  The  two  aspects  have  been  kept  separate  to  facilitate  the 
experimental  use  of  PROP,  but  it  will  be  normal  of  course,  if  the  user 
considers  ho  Li  near  tho  eccentricity  sinf^larity,  cay,  to  set  JELTYP  to  deal 
with  both  aspects  together. 

The  effect  on  the  computer  output  rust  be  mentioned.  For  the  line- 

printer  output  there  is  no  effect  at  all  on  the  final  results  (assu-ming  that 

final  results  could  be  obtained  if  the  two  aspects  of  the  singularities  were 

not  tackled'  );  in  particular  the  correlation  matrix  (S9,5)  still  gives  the 

correlations  between  standard  parameters.  For  the  card-punch  output,  tackling 

the  first  aspect  has  no  effect.  Tackling  the  second,  however,  leads  to  the 

punching  of  a  non-standard  covariance  matrix;  thus  if  partial  derivatives  with 

respect  to  fi  ,  w  -»  k  fi  and  M  -•  w  •*  k  n  have  been  introduced,  the 
o  o  o  o  o  o 

pxinched-card  covariamce  matrix  relates  to  these  parameters. 

6.f  Convergence 

The  aim  of  the  PROP  refinement  procedure  is  to  adjust  the  orbital  para¬ 
meters  in  such  a  way  as  to  minimize  the  sum  of  the  squares  of  the  weighted 
residuals,  a  residual  being  the  difference  between  am  observed  quamtity  amd 

the  corresponding  quamtity  calculated  from  the  dynamical  model;  symbolically 

2 

the  aim  is  to  minimize  w  R  .  where  R  -  ^  -  C  and  w.  R.  o,  c  denote 
respectively  weight,  residual,  observed  quantity  amd  calculated  quantity. 

The  procedure  is  am  iterative  one  and  a  test  for  convergence  is  made  at  the 
end  of  each  iteration. 
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The  convergence  test  is  based  on  the  value  of  a  quantity  c.  stored  as 
the  FORTRAN  variable  EPSLON.  This  is  the  estimated  standard  deviation  of  an 
observation  of  unit  weight;  if  the  total  number  of  observed  quantities  is  n 
and  the  number  of  orbital  parameters  under  refinement  is  N.  then 

c  -  {  r  w  /  (n  -  N)}^ 

there  being  n  -  N  df^grees  of  freedom. 

If  e  decreases  from  one  iteration  to  the  next,  the  refinement  pro¬ 
cedure  is  considered  to  be  converging.  When  the  decrease  is  less  than  one 
per  cent  the  convergence  is  deemed  to  be  complete,  so  long  as  the  process 
has  reached  the  zero  mode  (s6.^*)  sind  observations  are  not  still  being  rejected. 
If  the  dynamical  model  is  a  good  representation  of  the  actual  satellite  motion, 
and  if  the  observations  have  been  given  realistic  w'eights,  then  the  final 
value  of  E  should  lie  between  about  ''.7  and  1.5.  A  value  in  excess  of  ?  or 
h  should  be  viewed  with  some  disquiet. 

Due  to  the  non-linearity  of  the  system,  e  may  increase  instead  of 
decrease.  If  this  happens  twice  in  succession  the  refinement  procedure  is 
considered  to  be  divergent  and  is  discontinued.  However,  it  is  possible  for 
the  procedure  neither  to  converge  nor  to  diverge;  for  example,  c  may 
oscillate.  It  is  essential,  therefore,  that  the  PROP  user  set  a  maximum  to 
the  number  of  iterations  that  will  be  allowed.  This  is  done  by  one  of  the 
parameters  -  MAXITN  -  on  the  control  card, 

6.7  Rejection 

A  major  problem  in  the  processing  of  observational  data  in  any  field  is 
the  rejection  of  poor  data.  With  orbit  determination  the  basic  difficulty 
is  that  a  good  orbit  cannot  be  obtained  until  bad  observations  have  been 
eliminated,  while  they  cannot  be  judged  '*bad"  with  any  certainty  until  a  good 
orbit  is  available.  An  additional  difficulty  with  PROP  is  that  rejection  has 
had  to  be  made  entirely  automatic  -  it  cannot  be  left  to  the  user  to  decide 
at  the  end  of  one  iteration  how  to  reject  during  the  next  iteration. 

The  PR0P3  rejection  mechanism  has  the  following  three  main  features: 

(i)  its  basis  is  the  ratio  of  the  weighted  residual  associated  with 
an  observed  qusmtity  to  the  current  value  of  e  (where  an  observation 
consists  of  more  than  oi.e  observed  quantity  it  is  the  numerically  largest  of 
the  weighted  residuals  which  is  considered); 


(ii)  the  decision  as  to  which  observations  to  reject  is,  as  far  as 
possible,  postponed  until  the  end  of  the  Iteration; 

(lii)  rejection  of  an  observation  is  never  irrevocable  -  in  subsequent 
iterations  it  may  be  accepted  if  its  residuals  are  smaller. 

To  implement  featuws  (i)  and  (il)  the  observations  are  divided  into  four 
groups,  or  levels,  during  any  iteration.  Currently  (S5.2)  the  magnitudes  of 
the  weighted  residuals  for  all  observations  at  the  lowest  level  are  less  than 
y:-,  the  magnitude  of  the  largest  weighted  residual  for  each  observation  at 
the  second  level  is  between  y  and  at  the  third  level  it  is  between  ^ 
and  lOe,  and  at  the  highest  level  it  exceeds  ICP.  Cfbservatlons  at  each  level 
will  be  either  all  accepted  or  all  rejected.  The  lowest  level  is  always 
accepted  and  the  highest  always  rejected.  The  program  decides  at  the  end  of 
the  iteration  whether  to  accept  one  or  both  of  the  two  middle  levels;  the  basis 
for  this  will  not  be  described  here. 

Implementation  of  feature  (ili)  leads  to  the  slight  risk  of  an  oscillation, 
against  which  the  PROP  user  must  be  warned.  A  small  set  (say  two  or  three)  of 
observations  may  be  rejected  during  an  iteration,  and  the  effect  on  the  orbital 
parameters  and  c  may  be  such  that  on  the  next  iteration  they  have  to  be 
accepted,  after  which  they  are  rejected  again  etc. 

If  too  many  observations  have  been  rejected  at  any  stage  of  the  refinement 
process  it  becomes  pointless  to  continue.  The  PROP  user  should  specify,  by  the 
control  parameter  MINOBS.  the  minimum  number  of  observations  with  which  the 
process  will  be  allowed  to  continue  at  any  stage. 

6.8  The  ”previous-orbit"  facility 

It  sometimes  happens  that,  after  the  determination  of  a  set  of  orbital 
parameters  at  a  given  epoch,  further  observations  for  the  same  period  become 
available.  It  would  be  possible  to  make  a  fresh  determination,  using  the  old 
and  new  observations  together,  but  this  would  be  inefficient.  The  same  results 
could  be  obtained  by  using  the  new  observations  only.  If  the  PROP  user 

(a)  sets  the  control  parameter  NOTHER  to  a  non-zero  value  (as  explained 
below), 

(b)  takes  the  final  elements  of  the  previou.'^  run  as  initial  elements  for 
the  new  run,  and 

(c)  supplies  the  computer  with  "previous-orbit"  card.s,  which  give  the 
final  c,  the  number  of  degrees  of  freedom  and  the  covariance  matrix  of 
orbital  parameters  from  the  previous  run  (s8.5). 


If  N'^'IFP  t,‘ie  co*r}nil^r  vl  I !  n't.  r<*«rl  ;  r»v!  oue-ort t  rards.  If 
NOTHIP  1  1*.  r<*R!n  them,  taklnn  t*s«*  cr/arianc**  fr.alrlx  as  In  standard  forr. 

If  NOTJCF  ■  1  the  onputer  will  inter,  ret  the  covariance  matrix  as  having 
been  derived  from  a  non-standard  setting  of  JTI.TY!  In  the  jrevl'us  run 
(  St  .  and  r,“  .  2 ) . 

If  the  ;.revl ouB-orbl t  fa-l‘.:ty  is  use  1  MODE  should  always  be  set  to 
zero  since  the  initial  parameters  rust  be  "good".  If  .VODF  is  set  greater 
than  zero  PHOP  will  in  fact  Ignore  the  previous-orbit  data  until  MODE  has 
dropped  to  zero. 

There  is  one  situation  when  running  a  fresh  determination  with  all  the 
obsen’ations  may  lea!  to  slightly  different  results  from  using  the  previous- 
orbit  facility.  In  the  former  case  the  presence  of  new  observations  may 
cause  a  different  subset  of  the  old  obser/atl ons  to  be  rejected  in  comparison 
with  the  subset  previously  rejected.  In  the  latter  case  the  subset  of 
(implied)  rejections  cannot  var:/. 

7  PRELIMUiAFY  PPa’^ESSPlG  OF  OP"EFP>'ATIO?.‘: 

The  *'-0  main  stages  in  any  PPOP  run  are  the  input  of  obsen’atl ons  - 
subroutine  OBflf.’  -  and  the  differential  -'crre'-tlon  of  orbital  parameters  - 
subroutine  DL*''’3R.  In  between  these  is  a  subsidiary  stage  during  which 
certain  pr*’ 1  Iminar.'  '.'rocesslng  of  the  'bser.’ati  ons  Is  carried  out.  This 
subsidiary  stage  organizes  the  rotation  of  station  coordi nate.s  from  an 
earth-fixed  coordinate  system  to  the  standard  PROP  system.,  and  is  also 
responsible  for  the  -orrectlon  of  the  obser.-ntl  ons  for  refraction  llght- 
tlr.e  and  precession  and  nutation.  The  subroutine  responsible  for  all 
this  is  called  PROVES. 

7.1  Rotation  of  station  coordinates 

The  coordinates,  of  each  sensor  are  rotated  from  on  earth-fixed 
system  of  a.xes,  such  that  the  x-axls»  points  towards  the  Greenwich  meridian 
and  the  z-oxis  towards  the  north  pole,  to  t.-ie  standard  PRO!  system  (J*). 

The  rotation  is  about  the  z-axis,  \  an  amount  *  where  b  is  the 
"modified  sidereal  angle",  given  by 

7  1  '.  •  *'  «•  oQr,  T  \  '  ''2^2 .  '  : 


here  d  refers  to  'TTl  time  m.eas’urer}  in  VJT. 


2k 


2'i‘> 


No  »r -mini  '.a  Inker,  of  the  polar  notion  !,•.  of  the  motion  of  the  enrth’n 
rotation  ajila  relative  to  ita  cnjet.  The  reaultlnit  error  should  be  lerj  than 
15  *  In  itat'on  poaltlon  1.'-,  lees  than  '^".5  In  direction  of  the  local 
vertical . 

7.2  Refraction  correction 

.•llnce  all  methods  of  observation  depend  on  the  transjT  1  s.i  1  n  of  »-'ertro- 
mafcnetlc  radiation,  and  since  all  su^'h  radiation  Is  refracted  by  the  atmosphere. 
It  Is  necessary  to  make  appropriate  corrections  to  the  obser/atl ons .  For 
observation.*,  of  rl*:ht  ascension  and  declination  It  Is  not  appropriate  to  make 
the  total  f-.)rr''''t  1  or  from,  obser/ed  line  of  sight  to  true  line  of  sight  because 
tne  observations  ar**  male  relative  to  stars  which  are  themsc ive.*  seen  by 
refracted  radiation;  only  the  "paral  la'-tl  ^  correction"  Is  appropriate  In  this 
•'aae.  Whether  the  observed  quantities  were  orlKif'slly  arlruth  and  elevation  or 
right  aacenslon  and  declination.  It  Is  elevation  which  must  be  corrected  so 
that  In  the  case  of  original  right  ascension  and  declination  an  Initial  trons- 
foraatlon  Is  made. 

The  formulae  used  for  the  correction  of  elevation  are  the  sane  as  were 
used  In  the  old  Pegasus  progran^  but  are  repeated  In  the  next  two  paragraphs 
for  convenience.  The  correction,  whether  of  total  refraction  or  parallactic 
only  Is  a  function  of  the  observed  elevation  Itself  and  the  geocentric 
distance  of  the  satellite.  For  the  latter  quantity  a  value  Is  computed  from  the 
Initial  values  of  the  orbital  elements,  and  this  should  normally  be  sufficiently 
accurate . 

For  total  refraction  let  E'  be  the  observed  elevation  E  the 
corrected  elevation  and  let  r  and  F  be  the  geocentric  distances  to  the 
satellite  and  station  respectively.  Let  n  be  given  by 

tan  2n  ^  cot  E* 


and  1  »•  t 


5  5  7 

E'  -  (a  tan  n  *  b  ton'^  n  ••  c  ton'^  n  -•  d  tan  n  ) 


where  a  f'5^".515.  b  c  9  and  d  2^P".5n. 

[/•t  y  and  be  derived  from 
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r/> 


COB  1  “1.  /X)27^^5^  (B/r)  co«  E' 


and 


0 


'bo 


Then  E  la  Riven  by 


tan  E  '  cot  ^  -  (P/r)  cosec  ^ 

For  parallactic  refraction  ^  is  observed.  Let  be  Riven  by 

tan  cot  ^ 


and  let 


•  c  n 

kn  '  ^  ^  ‘  ^  ^  *  ''  tan^  n  ♦  d  tan  n 

where  a  -  95<''’.^2^  b  -  5^^".C95.  c  ^  175". 2;Q  and  d  -  °".?V  . 

Then  the  forrrulae  for  ■».  and  E  are  the  sorf  as  above. 

For  radio  obser/atlons  the  or-.ount  of  refraction.  E'  -  E.  is 
multiplied  by  a  factor  of  5.  (Thl.'.  factor  was  obtained  fror  Ref.®  and  is 
presu-nably  related  to  a  particular  fr**quency.) 

In  the  present  version  of  PROP  no  r'*fractlon  correction  is  nade  for 
observations  of  other  types.  In  parti'nilar.  ron^e  range  rate  and  direction 
cosines  are  not  corre^'ted.  It  is  be’leved  that  MlnitrncV.  observations 
(direction  cosines)  and  some  range  obsei-.-ations  are  corrected  for  refraction 
before  distribution. 

7. 5  Light-tlr.e  correction 

As  with  the  old  Pegasus  progra.-'.^  the  finite  speed  of  light  is  allowed 
for  by  correcting  the  reported  time  of  each  observation  back  to  the  time  at 
which  the  light  (or  radio)  waves  left  the  satellite.  This  requires  an 
approximate  value  of  the  topocentrlc  range  of  the  satellite  to  be  known.  If 
range  is  one  of  the  observed  quant Itle.s  it  is  used  d  rectly.  Otherwise,  the 
approximate  range  l.e  computed  from  the  initial  orbital  elements. 


26 


299 


7.^  Correction  of  right  aicenalon  md  declination  for  prtc«i«lon  and 
nutation 

Vfhether  right  eacenelon  and  declination  are  given  relative  to  the  mean 

equator  and  equinox  at  one  of  the  epochs  1^55-0.  .  1900.0  or  1950.0 

(the  epoch*  of  the  principal  atar  catalogues),  or  whether  they  are  given 

relative  to  the  equator  and  equinox  of  date,  a  trana  format  Ion  to  the  atandard 

PROP  coordinate  system  (S5)  la  carried  out.  Let  and  6^  be  the  quoted 

valuea  of  right  ascenalon  and  declination,  aaaumed  relative  to  a  atandard 

epoch,  let  o'  and  6*  be  their  valuea  relative  to  the  mean  equator  and 

equinox  of  date,  and  let  o  and  t  be  their  valuea  relative  to  the  stant'.ard 

PROP  ayatem.  The  transformation  from  o  and  t  to  a  and  6  is  made  in  three 

o  o 

atagea.  (Only  the  last  stage  la  Involved  If  the  anglea  are  given  relative  to 
the  true  equator  and  equinox  of  date.) 

Stage  1  The  first  transformation  is  from  and  to  o'  and  6'.  l,e. 

It  corrects  for  precession  between  the  date  of  the  atar  catalogue  and  the 
date  of  the  observation. 

Let  k  be  the  nuaiber  of  days  from  the  catalogue  epoch  to  195^.0. 
l.e.  5A698  27593  IB262  and  0  for  the  epochs  of  1855.0.  1075.0.  1900,0  and 
1950.0  respectively.  Let  d  be  the  days  from  1950.0  to  the  time  of  the 
observation  and  let  d*  =  k  *  d.  (For  an  observation  quoted  relative  to  the 
mean  equator  and  equinox  of  date  d'  0  and  k  *  -d.) 

Then  the  tranafomatlon  may  be  expressed  in  terms  of  the  standard^  ^ 
precesalonal  quantities  z  and  f’.  given  In  radians  by;- 

■=  (3.059557  X  lo""^  -  5.0781  X  10"^^  k)  d*  4  1.0986  X  10“^^  d»^  , 

z  -  K  *  2.B809  X  lo""*^  d'^  and 
o 

0  ^  (2,660^0  X  lo""^  4  5. ion  X  10"^^  k)  d'  -  1.5507  X  10"''^  d'^ 

The  required  transformation  is  given  by:- 


^COfl  ^  '  coo  <1  •' 
coo  f '  0 In  n • 

'  0 1  n  ft  '  / 
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a 

fCOB  Z  COB  *' 

-  sin  z 

-  COS  z 

s  i  n 

1  COB 

ft 

o 

cos  (a  *  K  )  \ 
o  o' 

sin  r  COB  '■ 

1 

COB  Z 

-  sin  z 

B 1  n 

0 

1 

CO.S 

1 

ft 

0 

sin  (a  *  K  ) 

0  o' 

\  sin 

cos 

♦ 

1 

1 

sin  ft  / 

o 


Stage  2  The  next  trano  format  Ion  Is  from  o'  and?'  to  values  relative  to 
the  true  equator  and  equinox.  These  nay  be  denoted  by  o"  and  ft,  since  the 
true  equator  is  the  final  PROP  equator  required.  This  transformation  corrects 
for  nutation . 

Let  and  be  arguments  of  the  moon's  orbit  such  that,  in 

radians. 


-  0.211  n  -  0.2422  X  10"^  d 

"hi  =  3.4955  ^  3.44  y  in""  d  and 

"^^'m  "  ■  2.2473  *  4.5004^1  x  Iq"''  d 

Then  the  principal  terns  of  A,  the  nutation  in  longitude  and 
the  nutation  in  the  obliquity  are  given  by:- 

Av  -  1  "I  ^  sin  n  -  0.(2  sin  2L,  O.lO  sin  20 

M  1*1  M 

-  ^.10  sin  2(F.,  ^w))  and 

M  M 

Ac  =  lo"^  (4.46  cos  4  0.27  cos  2L,  -  0. '•4  cos  20,, 

M  M  M 

-•  "^.04  cos  2(Fu  ^  “»))  ■ 

M  M 


The  required  transformation  is  given  by:- 
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/  coa  6  coa  o "  ' 

m 

I 

1 

-  Aa  coa  c 

-  A.a  Bln  c\ 

coa  ft  sin  Q* 

Aa 

coa  t 

1 

... 

\  aln  ft  ^ 

aln  t. 

Ax 

1  1 

where  iln  c  -  \59788l21  and  con  c  .  ). 91745694. 

Stage  5  The  final  tranifomatlon  la  from  a"  to  o .  This  arlaes  from  a 
rotation  within  the  true  equator  from  the  true  equinox  to  the  pseudo-equinox 
uaed  by  PROP, 

Define 


6.119  X  lo”*^  d  (in  radians) 
c> 


Then 


a  =  a"  -  41  -  A*  cos  c 

w> 


8  THE  IA1A  DECK 

8.1  Overall  description 

This  section  explains  how  to  prepsure  a  data  deck  of  punched  cards  for 
PR0P5.  l.e.  the  complete  set  of  cards  read  under  vairious  Input  formats  during 
the  various  runs  of  a  Job.  The  listing  of  a  typical  data  deck  is  attached, 
for  reference,  as  Fig. 4.  It  is  not  possible  to  give  instructions  for  the  pre¬ 
paration  of  a  complete  Job  deck  because  this  depends  on  the  job  description 
cards  required  by  the  particular  computer  being  used  whether  "source"  or 
"binary"  cards  are  being  reeui  for  the  program  itself,  whether  some  or  all  of 
the  program  is  held  on  magnetic  tape  etc.  (On  the  1907  computer,  if  the 
program  is  run  from  magnetic  tape,  as  overlaid  in  binary,  then  no  job  descrip¬ 
tion  cards  at  all  are  required.) 


Instructions  for  punching  will  not  pre-suppose  any  knowledge  of  FORTRAN. 
However,  the  actual  FORTRAN  input  formats  are  always  given  for  completeness. 

It  is  remarked  that,  according  to  A.S.A.  rules,  the  character  '*blank" 
(indicated  in  what  follows  by  the  letter  b)  is  completely  equivalent  to  zero. 
However  Atlas  is  known  to  object  to  embedded  blanks;  thus  it  is  better  to 
punch  680202  (for  the  date  I968  FEB  2,  read  as  a  6-diglt  integer)  than  to 
punch  68b2b2,  The  punching  of  plus  signs,  though  optional,  should  be 
eschewed,  since  "plus"  (unlike  "minus")  differs  as  between  the  Atlas/lBM 


I 


and  ICT  card  codes . 
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The  order  of  the  data  carde  for  a  nm  is 

(l)  control  card, 

(U'  epoch/identity  card  orbital  elements  (normally  on  five  cards). 

vlll)  previous -orb It  cards. 

(iv)  sensor  cards  blAnk  card. 

(v)  time-system  correction  cards. 

(vi)  observation  cards  -*  blank  card. 

It  will  be  seen  In  due  course  that  any  of  these  six  items,  except  the  firJt. 
may  be  omitted  In  appropriate  circumstamces .  The  blank  card  in  item  (iv)  is 
actually  read  as  a  sensor  card  and  indicates  tl-,at  no  more  sensor  cards  follow. 
With  Item  (vi)  the  situation  is  more  complicated.  A  'set*  of  observation 
cards  is  made  up  of  one  or  more  subsets,  as  explained  in  section  8,8;  eacn 
subset  ends  with  a  blank  card  and  so  item  (vi)  ends  with  two  blank  cards.  Thus 
a  run  ends  in  two  blank  cards  whenever  item  (vi)  is  present  (but  in  no  blank 
cards  if  both  items  (iv)  and  (vi)  are  absent). 

The  complete  data  deck  for  a  Job  consists  of  the  data  cards  for  the 
various  runs,  followed  by  a  blank  card;  this  is  read  as  a  control  card  and 
indicates  that  there  are  no  further  runs.  Assuming  item  (vi)  is  included  in 
the  last  run.  a  complete  data  deck  must  end  with  three  consecutive  blank  cairds. 

A  general  warning  on  punching  must  be  made.  Integers  must  always  be 
punched  as  far  right  as  possible  in  their  respective  fields,  since  they  will 
otherwise  be  misinterpreted.  A  mixed  number,  on  the  other  hand,  may  be 
punched  anywhere  in  its  field,  provided  that  a  decimal  point  is  included. 

8.2  The  control  card 

The  control  card  for  each  run  contains  twelve  parameters  (all  integer), 
which  control  the  working  of  the  program  during  that  run.  It  also  contains 
optional  re-settings  of  four  of  the  geophysical  and  program  constants  (one 
integer  and  three  real).  The  input  format  is  (15I5'  which,  since 

every  field  is  five  characters  in  w’dth  makes  it  easy  to  punch.  The  inter¬ 
pretation  of  the  twelve  control  parameters  and  four  optional  constants  is 
given  below,  with  some  indication  of  the  values  that  should  normally  be 
punched.  The  parameters  and  constants  are  referred  to  by  their  FORTRAU- 
variable  names  as  used  in  the  program.  One  possible  punching,  for  reference 
is  as  follows  : 
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( - 1 - r  “  ■' - ’ - r  - - T  *  ■  ^ - - - 1 - - - : - - - 1 - -  i 

222241  m  4bbbb0bbbb1  bbbb5bbbb7bbb20bbb1 1  bbbb2bbbb0bbbb2bbbb1  bbbb9bb2 . 1  bbbbbbbbbb 

nrrrr  (colJ.1-5)  This  Indicates  which  elements  are  present  In  the  orbital 

model.  The  five  digits  of  KKKXK  are  associated.  In  order, 
with  the  polynomials  which  represent  e.  T.  tj  and  m;  for  each  digit  K. 
the  degree  of  the  corresponding  polynomial  Is  K  -  1.  E.g.  If  KKKKX  ^  22224, 

then  e,  T,  are  linear  and  ft  Is  cubic;  the  elements  present  in  this 

case  would  be  e^,  .  1^.  Lj.  0^.  0^.  M^.  Mj.  The  standard 

model  (S4.2)  gives  a  KKKKK  of  22223  or  22224,  (N.B.  If  account  is  taken 
of  the  first  point  discussed  In  section  4.ii  it  Is  logical  to  make  the  degrees 
of  the  e.  1.  n  and  polynomials  one  less  than  the  degree  of  the  M  poly¬ 
nomial.  so  that  the  standard  values  of  KKKKK  become  22223  and  33334,  The 
difference,  in  orbit  determination  over  a  few  days.  Is  negligible.) 

MUflM  (cols.6-10)  This  Indicates,  similarly,  the  elements  which  are  (float¬ 
ing)  parameters  to  be  refined.  E.g,  if  MMMMM  =  11114 

then  e  .  i  ,  n  ,  ,  M  .  M,,  Mr,  WJd  M,  are  parameters;  if  MMMMM  =  10l02 

0000012  \5 

then  only  e^,  M^  and  M^  are  parameters.  The  standard  model  gives  an 
M^#<MM  of  11113  or  11114.  corresponding  to  the  standard  values  of  KKKKK. 

Clearly  no  m.  where  m  is  a  digit  of  MMMMM.  may  exceed  the  corresponding 
K.  and  at  least  one  m  must  be  non-zero. 

MODE  (cols. 11-15)  This  indicates  the  mode  in  which  refinement  is  to  take 

place  on  the  first  iteration  of  the  differential-correction 
process.  If  the  initial  estimates  of  the  orbital  parameters  are  believed  to 
be  good,  punch  MODE  =  0;  if  they  are  uncer  .n  in  quality,  punch  1;  if  they 
are  believed  to  be  rather  poor,  punch  2.  etc,  (S6.4). 

NEWS.4T  (cols. 16-20)  This  indicates  the  way  in  which  the  epoch  and  initial 

elements  for  the  run  are  provided.  There  are  five 
possibilities,  as  follows  :- 

IIEWSAT  =  0:  epoch  and  elements  will  be  taken  directly  from  the  six  cards  which 
immediately  follow  the  control  card;  l.e.  epoch/identity  and  element  cards 
cannot  be  omitted. 

NEWSAT  =  1  :  epoch/identity  and  the  five  element  cards  are  present,  as  in  the 
pirevious  case,  but  the  run  is  actually  to  be  associated  with  a  different  epoch 
(usually  a  few  days  later)  and  the  input  elements  are  to  be  used  to  predict 
initial  elements  at  the  new  epoch. 
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NEWSAT  ‘=2:  no  epoch/identity  or  element  cards  are  provided;  initial 
elements  for  the  run  are  to  be  obtained  by  prediction  from  elements  stored 
in  the  computer  at  the  end  of  the  run  Juat  completed. 

NEWSAT  5:  this  is  the  same  as  the  previous  case  except  that  a  special 
card  is  provided  from  which  new  values  of  the  exclusive  elements  e.j ,  i^ , 
and  I  are  taken. 

NEWSAT  this  is  similar  to  the  first  case  (NEWSAT  ),  but  after  the 

epoch/identity  card  appears  a  single  card,  containing  the  six  components  of 
initial  position  and  velocity,  instead  of  the  five  element  cards;  these  will 
at  once  be  converted  to  mean  orbital  elements  with  any  additional  (ex'^lusive) 
elements  in  the  model  set  to  zero. 

(N.B.  (1)  If  NEWSAT  -  2  or  5  there  must  be  at  least  one  preceding  run  in  the 
data  deck,  and  this  must  be  for  the  same  satellite  as  the  new  run,  since  other¬ 
wise  prediction  from  the  preceding  run  would  have  no  meaning.  If  NEWSAT  , 

1  or  thcr^  may  or  may  not  be  preceding  runs  and,  if  ther'--  are,  the  new  run 
may  be  for  the  same  satellite  or  a  "new  satellite". 

(2)  There  is  no  NEWSAT  value  \hich  covers  position  and  velocity  plus 
prediction . ) 

MJDINC  (cols. 21 -25)  This  is  used  if  prediction  is  required  (NEWSAT  ^1,2 

or  5).  It  gives  the  increment  in  Modified  Julian  Days, 
to  be  applied  to  the  current  epoch  (MJDOCH)  to  bring  it  to  the  epoch 
required  for  the  new  run.  E.g.  to  predict  forward  three  days  punch  to 
predict  back  a  week  punch  -7. 

MAXITH  (cols. 26-50  This  specifies  the  maximum  number  of  iterations  that 

will  be  allowed  by  the  refinement  process.  E.g.  to 
stop  a  particular  run  after  ten  iterations,  whether  the  process  has  converged 
or  not,  punch  lo.  (if  the  PROP  rcflr.tment  process  is  going  to  converge  it 
will  rarely  require  more  than  ten  iterations.) 

MINOBS  (cols. 31-55)  This  specifies  the  minimum  number  of  obser/ations  with 

which  the  refinement  process  will  be  allowed  to  continue 
at  any  stage.  Thus  it  might  be  thought  intolerable  to  have  miore  tnan  a  tnird 
of  the  observations  rejected  at  any  stage,  in  which  case  with  5  obser/aflons 
initially,  say,  MINOBS  would  be  21, 

JELTYP  (cols.3t-^  )  This  indicates  whether  the  standard  elements 

e  ,  i  ,  n  ,  ,  M  are  used  throughout  the  run.  or 

o  o  o  o  o 

whether  non-singular  elements  are  introduced  at  some  point.  JELTYP  is  best 
thought  of  as  a  two-digit  number,  of  which  the  first,  or  tens  digit  (wnicn 
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app«art  In  coluam  59)  Bay  be  dtnoted  by  JELl  and  the  second,  or  units  digit 
(which  appears  in  column  4o)  nay  be  denoted  by  JTLS.  J^il  and  JEl£  are 
treated  independently  by  the  program  and  since  each  can  have  the  value  0.  1. 

2  or  5  there  are  16  possibilities  for  JELTYP.  JTLI  is  associated  with  the 
transforeation  of  parameters  to  avoid  ill-conditioned  matrices  (S6.3)  as 
follows 

JEL1  ■  0  :  no  transformation; 

JELl  =  1  :  transformation  to  •'  and  M  +  '  ,  when  e  is  close  to  zero; 

o  o  o  o 

JELl  =  2  :  transformation  to  n  and  +  k  0  .  when  i  is  close  to  or 

o  o  o  0 

18T; 

JELl  =  3  !  tr^jisformatlon  to  0  ,  m  +  k  0  and  M  +  ’■’  +  k  0  .  when  both 

o  o  o  o  0  o 

e  and  i  are  near  to  singular  values, 
o  o 

JEL2  is  associated  with  the  temporary  transformation  of  elements  to  improve 
the  convergence  of  the  PROP  refinement  process  (S6,5)  as  follows 
JEl£  =  0  :  no  transformation; 

JEI^  =  1  :  transformation  to  e  cos  e  sin  m  and  M  +  ■'  .  when  e  is 

oooooo  0  • 

close  to  zero; 

JEL2  =  2  :  transformation  to  sin  i  cos  0  and  sin  i  sin  Q  .  when  i  is 

o  O  0  0  0 

close  to  0*  or  l8o*; 

JEI2  =  5  ’  transformation  appropriate  when  both  e^  and  i^  are  near  to  singular 
values . 

Experience  shows  that  the  choice  of  a  suitable  JEL2  is  more  important  than 
JELl;  however,  a  simple  procedure  is  to  choose  the  same  value  for  both,  so 
that  JELTYP  is  o,  Tl.  22  or  55. 

ROTHER  (cols. 41 -45)  This  indicates  whether  or  not  previous-orbit  cards  are 
Included  in  the  data  deck  for  the  run: 

ROTHER  =  ')  :  previous-orbit  cards  are  not  included; 

ROTHER  =  1  :  previous-orbit  cards  are  included  and  JELl  =  0  for  both  the  previous 

and  new  iruns ; 

ROTHER  -  2  :  previous-orbit  cards  are  included  euid  JELl  ^  1  for  both  the  previous 
and  new  runs; 

ROTHER  -  3  :  previous-orbit  cards  are  included  and  JELl  =  2  for  both  the  previous 
and  new  runs; 

ROTHER  -  4  :  previous-orbit  cards  are  included  and  JELl  =  3  for  both  the 
previous  and  new  runs. 

R.B.  previous-orbit  cards  may  only  be  included  if  REWSAT  =  0. 
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ITIMEC  (cola. 46-50)  This  indicates  what  time-system  correction  cards  if 

any,  are  included  in  the  deck. 

ITIMEC  =  0  :  no  time-correction  cards  are  Included  for  this  run; 

ITIMEC  =1  :  a  single  WWV  time-correction  card  is  included; 

ITIMEC  =  2  :  a  pair  of  VWV  time-correction  cards  is  included; 

ITIMEC  =  10  :  an  A1  time-correction  card  is  included; 

ITIMEC  =  n  :  an  A1  and  one  WWV  card  are  included; 

ITIMEC  =12  :  an  A1  and  two  WWV  cards  are  included. 

ISENSR  (cols. 51 -55)  This  relates  to  the  reading  and  printing  of  sensor-caard 

data. 

ISENSR  =  0  :  sensor  cards  are  present,  but  information  is  not  to  be  printed 
after  the  reading  of  each  card. 

ISENSR  =  1  :  sensor  cards  are  present,  and  information  about  each  card  i^ 
to  be  printed  after  reading. 

ISENSR  =  2  :  sensor  cards  are  not  present,  and  the  subroutine  for  reading 
sensor  cards  is  to  be  by-passed  in  the  present  run;  PROP  will  use  sensor  data 
provided  by  an  earlier  run  in  the  same  Job.  (This  option  is  quite  independent 
of  the  NEWSAT  option;  so  long  as  the  MAXSTA  limit  is  not  exceeded  (S5.2)  it 
is  possible  to  pat  all  sensor-cards  for  a  job  into  the  first  irun  and  set 
ISENSR  =  2  for  all  the  other  runs . ) 

lOBSNS  (cols. 56-60)  This  relates  to  the  reading  and  printing  of  data  from 

observation  cards. 

I'OBSNS  =  0  :  observation  cards  are  present,  but  information  is  not  to  be 
printed  after  the  reading  of  each  card; 

lOBSNS  =  1  :  observation  cards  are  present,  and  information  about  each  card 
is  to  be  printed  after  reading. 

lOBSNS  =  2  :  observation  cards  are  not  present  and  the  subroutine  for  reading 
them  is  to  be  by-passed  in  the  present  run;  PROP  will  use  the  same  observations 
as  for  the  preceding  run  in  the  same  job;  preliminary  processing  will  be 
restricted  to  resetting  the  initial  rejection  levels.  (This  option  will 
rarely  be  selected,  but  occasionally  it  is  desirable  to  run  the  same 
observational  data  twice  in  succession  with,  say.  two  different  orbital 
models.) 

lOBSNS  is  the  last  control  parameter  on  the  control  card.  The  remaining 
columns  are  used  for  such  of  the  constants,  described  below,  as  it  is 
desired  to  reset.  A  blank  (or  zero)  field  for  any  of  these  constants  means 
that  the  "standard  value"  will  be  used. 


L  (cols. 61-65)  This  is  the  index  of  the  last  zonal  hannonic  to  be  used  bj’ 

PROP;  i.e.  J„,  J,,  ....  J  will  be  used.  The  standard 
^  j  L 

value  is  15  ®Jid  the  maximum  permissible  value,  unless  the  projntun  is  recom¬ 
piled  (S5.1)  is  L  =  16.  A  new  specification  will  apply  to  the  current  run 

only  and  PROP  will  simply  use  the  BLOCK  DATA  J's  up  to  J  . 

L 

EJ22  (cols. 66-70)  This  is  10  Jg  g  where  2  standard  tesseral 

harmonic.  The  present  standard  value  Is  1.^.  Tf  it  is 

required  to  be  zero  a  very  small  non-zero  value  must  be  punched.  Fairing 

most  of  a  PROP  run  EJ22  holds  the  value  of  ^  itself,  instead  of 

6 

10  jj  j. 

DEKSCH  (cols. 71 -75)  This  designates  H.  the  value  of  the  density  scale 

height  of  the  earth*s  atmosphere  at  a  height  '-H 
above  perigee,  in  km.  The  PR0P3  standard  value  is  25  (km)  and  alternative 
values  are  available  from  Fig.5>  taken  from  Ref.n. 

EPSLON  (co1s.76-Bi)  This  specifies  the  initial  value  of  r  (06.6),  the 

convergence  parameter.  The  standard  value  is  2"io.o. 

The  three  most  likely  x*easons  for  setting  a  different  value  are:- 

(ij  The  observations  to  be  used  in  the  run  are  of  very  hiph  nominal  accuracy 
(e.g.  Baker-Nunn  observations)  and  hence  their  weighted  residuals  will 
be  numerically  large  in  comparison  with  their  unweighted  residuals.  If 
the  initial  values  of  the  orbital  parameters  are  poor  the  initial  c 
must  be  much  larger  than  2'^^  or  else  all  the  observations  may  be 
rejected. 

(ii)  The  initial  values  of  the  orbital  parameters  are  thought  to  be  so  good 
(because  they  have  been  obtained  from  prediction,  for  example)  that  it 
would  be  unfortunate  if  an  observation  with  a  weighted  residual  of  50, 
say,  were  accepted  on  the  first  iteration.  To  ensure  rejection  of  such 
observations  a  much  smaller  initial  c  is  required.  The  obvious 
application  of  this  is  to  prevent  wasted  iterations  when  there  is, 
say,  a  single  "rogue"  observation  which  will  be  re.iected  at  the  end 
anyway,  but  there  is  another  application  which  should  be  mentioned.  If 
there  is  a  large  number,  more  than  55^  say,  of  poor  observations  the 
rejection  process  may  fail  completely  and  the  run  may  end  in 
"convergence"  to  a.i  unacceptable  final  c  with  none  of  the  poor 
observations  rejected.  This  cam  sometimes  be  avoided  by  starting  with 
an  c  as  low  as  2.  \  say.  Then  all  or  most  of  the  poor  observations 
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are  Imnedlately  reacted  at  the  hl^heit  level  (S6,7);  any  obaervatloni 
that  Btar'  \y  belnf?  wronp.ly  re.iected  will  be  accepted  In  later 
Iterations  if  all  goeo  well. 

(iii)  Tt  Ir  desired  to  effect  convergence  after  only  one  Iteration.  Ai  an 

exa-Tj^le  of  when  this  may  occur,  suppose  that  a  perfectly  succeaiful  run 
has  be»‘n  carried  out,  ending  with,  say,  EPSLON  ■  1,231,  It  li  then 
realised  that  due  to  uae  of  the  wrong  value  of  JELTTT  the  output 
covariance  cards  are  worthless.  To  repeat  the  run.  with  minimum  waate 
of  computer  time.  It  is  onb'  necessary  to  start  with  the  final  elements 
,'ust  obtained  and  an  initial  EPSLON  of,  say,  1,235.  since  the  convergence 
criterion  will  then  be  satisfied  (S6.^)  after  a  single  iteration, 

(N.n,  Tt  would  be  necessary  in  this  case  to  omit  all  the  observations 
which  were  rejected  by  the  final  Iteration  of  the  earlier  run.  since 
otherwise  their  fresh  rejection  would  prevent  single-iteration 
convergence . ) 

Tt  is  probably  sound  practice,  when  analysing  the  orbit  of  the  same  satellite 
at  regular  intervals,  to  start  all  runs  with  aji  EPSLON  of  5,0.  10,o  or  20.0. 

n. 5  The  epoch/identity  card 

An  epoch/identity  card  is  required  when  the  control  parameter  NEWSAT 
has  the  value  0.  1  or  ^ .  Tt  contains  three  items:  the  epoch,  stored  as  the 
(integer)  variable  MJDOCH;  the  satellite  identification  number,  stored  as 
the  (integer)  variable  IDElfT;  and  the  satellite  name,  stored  as  the 
(integer-Hollerith)  array  NAME(^).  The  input  format  Is  (l6,  5X,  17.  ^A4) 
and  a  typical  punching  for  reference  is  'b39^911>bb670^2cf-bAPIELb^bbbbbb’~ 
etc,  b, 

MJDOCH  (cols. 1-6)  This  specifies  the  Modified  Julian  Day  niunber  of  the 

epoch  at  which  the  orbital  elements  of  the  satellite  - 
which  will  be  punched  on  the  five  cards  (or  -ngj.e  card  if  NEWSAT  =  A)  follow¬ 
ing  the  cpoch/ldentity  card  -  are  defined.  If  NEWSAT  =  this  is  the  same  as 
the  epoch  with  which  the  run  is  to  be  associated.  MJDOCH  may  be  punched  as 
the  MJD  number  itself.  e,g,  b39891  for  19^8  FEB  5-  or  as  a  year-month-day 
number.  68'2")5,  (The  computer  interprets  the  number  read  by  checking  whether 
it  is  smaller  or  larger  than  I'looon.) 

IDENT  (cols  .1 'i-l6)  This  is  the  7-diglt  designation  used  by  the  Smithsonian 

Astrophyslcal  Observatory,  It  is  used  by  PROP  solely 


for  checking;  thnl  -AO  Baker-^J’inn  (and  I'ertaln  othei  )  oh n *•  rvn ♦  1  on. '•>  ar**  for  th<* 
right  Batelllt^,  [Vnotlng  the  7  diglta  by  YYT.LIiNTl  YY  given  the  year  of 
launch,  LLL  the  number  of  the  launch  within  that  year,  and  h7)  the  number  of 
the  object  within  that  launch,  F.g.  the  satellite  Ariel  ^  (1*'»7-^*2Al  Is 
dealgnated  ^7014201  and  the  rocket  of  that  satellite  h<^7-'t?R^  1 -c  ^70^202, 

NAME(^ )  (cols.  17-52)  This  may  be  punched  arbitrarily,  or  omitted.  Tt.s  only 

funetlon  Is  one  of  reference,  and  the  1^  allocated 
columns  are  simply  copied  to  the  printed  output,  A  natural  punching  would 
be  for  example.  -bARIELb^bbbbbbb  or  -bARTEI^b 5b ROCKET. 

Orbital  elements 

9.^.1  Five-card  orbital  elements* 

These  cards,  which  must  always  be  preceded  by  an  epoch/identity  card, 
are  associated,  one  each  In  order,  with  the  basic  elements  e,  i,  and  M. 

They  are  read  by  the  subroutine  ELRF.AD.  The  read  format  is  (IS,  IX,  T1, 
lx,  2F12.0,  itFn2.i|)  and  this  is  fully  compatible  with  the  lun^hed-card  out¬ 
put  of  (refln^'d)  elements  which  takes  place  at  the  end  of  a  PROP  ran.  This 
fact  Is  very  useful  when,  for  example,  orbital  elements  for  a  given  satellite 
are  being  computed  at,  say,  three-day  intervals,  with  .iobs  prepared  from 
batches  of  half-a-dozen  runs;  element  cards  which  have  beer  output  at  the  end 
of  the  last  run  of  one  Job  may  be  put  straight  into  the  data  deck  for  the 
first  run  of  the  next  ,iob,  with  NEWSAT  -  1  and  MTDITIC  -  5  on  the  new  control 
card  of  course. 

In  the  explanation  below  it  is  convenient  to  take  M  as  a  topical 
element  and  to  describe  the  punching  of  the  associated  card.  A  typical  punch¬ 
ing  is 


3qPrc|/5bbbb-n5.^2'4l4bbbSl<51  .0^^'^%bbbbb0 .0'4Q2etc(blank) 


i'ols.1-5  contain  MJDOu'H,  which  denotes  the  same  epoch  as  on  the  epoch/ 
identity  card  (n^.5)  and  must  be  punched  in  the  MJD  form  (or  else  omitted  - 
see  below);  col.f  contains  a  solidus  (optional);  col. 7  contains  N.  the  serial 
number  (1  to  5)  of  the  element  card,  e.g.  5  in  the  case  of  the  M  card,  (or 
is  omitted  if  MJDOCH  is);  and  col.^  is  always  blank.  These  first  eight  columns 
permit  a  program  check  that  the  element  cards  are  in  the  right  order  and  that 
they  correspond  to  the  epoch  given  by  the  epoch/identity  card.  If  the  check 

"riic.;o  nr<-  '!U'tp  different  fr.-M  the  five-r-ard  elements  issued  by  the 
U.ll.A.F.  in  their  Ppti -ctrack  bulletins. 


falla,the  run  -  and  indeed  the  real  of  the  job  -  la  abandoned.  The  check 
will  be  by-paaaed  for  any  card  for  which  cola, 1-8  are  blamk;  the  PROP  uaer 
may  take  advantage  of  thla  fact,  but  thla  la  not  recommended. 


Cola. 9  onwarda  of  the  M  card  contain,  in  flelda  of  width  12,  the 
coefflclentn  In  the  polynomial  (of  degree  -  1 )  which  repreaenta  M 
In  the  orbital  model.  Here  is  given  by  the  fifth  digit  of  the  control 

parameter  KKKKK.  E.g.  if  KKKKK  -  22225,  then  =  5.  *md  if  the 
quadratic  polynomial  repreaenting  M  Is 


-115.6241J  4  5i45l.o'4^8  t  4  0.0^92  t‘ 


the  M  card  could  be  punched  as  In  the  typical  example  above.  Here  t  is 
time  in  day.s  from  epoch  and  M  is  In  degrees. 

Punching  of  the  e,  i,  H  and  '  cards  similar,  but  it  must  be 

remembered  >  hat  the  dominant  term.s  of  e.^ .  i.^,  9.^  and  must  not  be  included 
i.e.  the  punched  values  are  of  the  exclusive  elements  (S^.5).  These 
exclusive  elements  will  often  be  zero,  in  which  case  cols. 21-52  may  be  left 
blank,  since  a  blank  field  Is  read  as  0.0. 

Units  are  degrees  (for  angles)  and  days  (for  time);  e.g.  e  is  non- 

2 

dimensional.  i.^  is  in  deg/d  and  9^  iri  deg/d  . 

One-card  position  and  velocity 

This  card  provides  an  alternative  to  the  input  of  five-card  elements 
and  is  used  when  NEWfAT  It  must  always  be  preceded  by  an  epoch/ 

identity  card. 

The  main  quantities  on  the  card  are  the  three  coordinates  of  epoch 
position  (x,  y,  z )  in  km  and  the  three  coordinates  of  velocity  (x,  y,  z) 
ill  km/sec.  The  read  format  is  (l*^.  ^X,  ^FlP.c);  cols. 1-5  provide  an 
optional  check  of  MJDOCH,  as  in  section  “^.^.l,  and  cols. 9-20,  21-52  etc,  are 
used  for  x,  y  etc. 

The  card  is  read  by  the  subroutine  PVTREAD,  which  converts  the 

coordinates  of  position  and  velocity  into  the  standard  mean  elements  e  , 

o 

i  ,  9  ,  ,  M  and  M,  . 

o  o  o  o  1 

8.1|.5  One-card  orbital  elements 

This  card  is  required  in  the  special  situation  which  arises  when 
NEWSAT  -  5.  There  will  be  no  epoch/identity  card  preceding  it. 
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It  is  recalled  that  NEWSAT  -  3,  like  NEWSAT  -  2.  Inv'^lves  the  pre¬ 
diction  of  elements  from  values  held  in  the  computer  before  the  beginning 
of  the  nin .  If  non-zero  values  of  the  exclusive  elements  e^.  i^,  ajid 
computed  outside  PROP,  are  being  supplied  for  each  run.  these  have  to  be  read 
in  to  supersede  predicted  values.  This  is  the  purpose  of  the  NEWSAT  -  5 
option.  It  is  emphasised  that  an  epoch/identity  card  is  not  required  for  this 
case. 

The  one-card  elements  are  simply  the  exclusive  values  of  e^  i.^, 
and  .  The  format  is  (4F20,0)  and  the  units  are  degrees  (for  angles)  and 
days  (for  time). 

RESDEL(1  )  (cols.1-2n)  This  is  the  exclusive  e.^ , 

RESDEL(2)  (cols. 21 -4n)  This  is  the  exclusive  i.^ , 

RESDEL(3)  (cols. 41 -6t)  This  is  the  exclusive  0.^. 

RESDEL(4)  (co1s.61-8o)  This  is  the  exclusive  i^>.j , 

8,5  Previous-orbit  cards 

The  significance  of  these  cards  has  been  explained  in  section  6.3.  They 
consist  of  a  single  card,  described  in  the  next  paragraph,  plus  a  set  of  cards 
which  contain  the  elements  of  the  covariance  matrix  as  punched  at  the  end  of 
the  previous  run,  (N.B.  'Previous  run*  has  the  meaning  implied  by  section 
6.0  and  does  not  relate  to  the  preceding  run.  if  any,  in  the  deck  of  cards 
for  the  current  Job.)  The  orbital  model  must  be  the  same  as  in  the  previous 
run  and  contains,  say,  I  parameters  given  by  summing  the  digits  of  the 
control  parameter  MMMMM,  The  covariance  cards  are  read  under  format  (8E10.5); 
if  f  8  there  will  be  I  such  cards,  if  9  f  16  there  will  be  21 
such  cards,  and  if  17  ^  f  20  (its  msiximum  possible  value)  there  will  be 
5i  cards , 

The  card  which  precedes  the  covariance  cards  contains  two  quantities, 
read  under  the  format  (l5,  Fl5.0). 

NOODOF  (cols  .1-5)  This  if  the  number  of  degrees  of  freedom  for  the  last 

iteration  of  the  "previous  run",  given  by  the  printed  out¬ 
put  for  that  run. 

OEPSLN  (cols. 6-20)  This  is  the  final  value  of  c  (s6.6)  for  the  previous 


run. 
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0.6  Senaor  cards 

Up  to  50  sensor  (observlng-station)  cards  may  be  Included  in  the  data 
deck  for  a  run.  (This  upper  limit  may  be  changed  if  desired;  see  section 
5.2.)  They  ore  read  by  the  subroutine  FENSIN,  assuming  the  control  para¬ 
meter  ISEINSR  is  less  than  2.  The  sensor  cards  must  be  followed  in  the  deck 
by  a  blank  card;  PROP  exits  from  SENSIN  as  soon  as  this  is  encountered. 

A  sensor  card  contains  information  about  the  sensor*s  location  (in 
either  geographical  or  Cartesian  coordinates),  its  reference  niunber  and  its 
name.  The  format  is  (ll,  19,  2X,  ^A4,  F12.0,  4F10.0)  and  typical  examples 
of  sensor  cards  are 

C22) 

Sbbbbbbl^^INkFIELDbbbbbb^. . b'bb5905 .15  3^-48 . 4o4\)b4964 . 71 1 


'bbbbbbi  000bbABAKANbbbbbbbbbbbbb9l  . 45855b'bb55 .721 67’bbb247 . Obb'b .  .^  . .  l  .b 


The  interpretation  of  the  quantities  to  be  punched  on  a  sensor  card  is  given 
below. 

MARKER  (col.1 )  This  specifies  the  type  of  information  provided  by  the  rest 
of  the  card.  There  are  four  possibilities 

MARKER  =  0:  geographical  coordinates  (longitude,  latitude  and  height)  are 
supplied  amd  they  refer  to  the  Fischer  spheroid^ 

MARKER  =  1  :  geographical  coordinates  are  supplied  and  information  which 
defines  the  reference  spheroid  is  also  supplied. 

MARKER  =  2:  Cartesiaui  coordinates  are  supplied  and  they  refer  to  the 
Fischer  spheroid. 

MARKER  =  5:  both  geographical  coordinates  (longitude  and  latitude  only) 
and  Cartesian  coordinates  are  supplied.  (This  effectively  means  that  the 
sensor  need  not  be  located  anywhere  near  the  surface  of  the  earth.) 

NOSTAT  (cols. 2-10)  This  is  the  station  number.  It  should  adhere  to  a 

recognized  numbering  system  such  as  the  COSPAR  system. 
E.g.  for  Villa  Dolores  punch  bbbbb90n  . 

NAME  ST  (cols  .15-28)  This  is  the  station  name,  stored  as  an  away  of  fo\ir 

Hollerith  elements.  Its  only  function  is  to  appear 
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with  the  output  printing  (if  ISENSR  -  1)  for  reference  purposes.  Punch,  for 
example,  VILLAbDOLORESbbb. 

3TAL0N  (  cols. 29-^0)  If  MWyCER  ♦  2  thin  is  the  longitude  in  degrees,  positive 

to  the  east  of  Greenwich.  If  MARKER  =  2  it  is  ignored. 

3TALAT  (  cols. ^1-50)  If  MARKERS  2  this  is  the  latitude  ir.  degrees,  positive 

for  the  northern  hemisphere.  If  MARKER  =  2  it  is 

ignored. 

STAHGT  (cols. 51 -6o)  If  MARKER  =  0  or  1  this  is  the  height  in  metres  above 

the  reference  .spheroid.  If  MARKER  -  2  or  5.  it  is  X 
In  km,  where  X  is  the  station ‘n  coordinate  measured  from  the  centre  of  the 
eaurth  towards  the  Greenwich  meridian. 

AEARTH  (cols. 61 -70)  If  MARKER  =  0  this  is  ignored.  If  MARKER  -  1  it  gives 

the  semi-major  axis,  in  km,  of  the  required  reference 
spheroid,  e.g.  637Q.38Q  for  the  Hayford  spheroid.  If  MARKER  -  2  or  3  it 
gives  Y,  where  Y  is  measured  towards  the  9  °E  meridian. 

HEARTH  (cols  .71 -80)  If  MARKER  =  0  this  is  ignored.  If  MARKER  =  1  it  gives 

the  semi -minor  axis  of  the  reference  spheroid,  e.g. 
6358.912  for  the  Hayford  spheroid.  If  MARKER  =  2  or  3  it  gives  Z,  measured 
towards  the  north  pole. 

8.7  Time-system  correction  cards 

Up  to  three  cards  for  correcting  the  time  system  (S3)  may  be  called  for 
by  the  program,  depending  on  the  value  of  the  control  parameter  ITIMEC  (S8.2). 

If  an  Al  (atomic  time)  correction  card  is  required,  it  must  come  first. 
It  contains  a  single  quantity,  as  follows,  reeui  under  format  (FlO.O). 

AlTC  (cols.l-lO)  This  gives  the  correction  to  Al  time  to  give  UTl ;  i.e.  it 

is  the  epoch  value  of  OTl-Al.  It  should  be  punched  as  a 
time  in  seconds,  with  decimal  point  included. 

If  a  single  WWV  correction  card  is  required,  this  contains  a  single 
quemtity,  read  under  format  (FlO.O),  as  follows. 

WWVTC  (cols  .1-10)  This  gives  the  correction  to  WWV  time  to  give  UTl.  Its 

punching  is  the  same  as  for  AlTC. 

If  a  pair  of  WWV  cards  is  required,  the  first  simply  contains  WWVTC  as 
above.  The  second  card  defines  a  discontinuity  in  the  WWV  system  and  contains 
two  quantities,  read  under  format  (16,  FlO.O),  as  follows. 
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JMPMY  (cols .1-6)  This  gives  the  date  of  the  discontinuity,  Oh  understood. 

Like  an  epoch  (sS.J)  it  may  be  punched  in  two  possible 

ways. 

VfWVTC2  (cols. 7-1 6)  This  gives  the  correction  to  be  applied  after  the 

discontinuity,  the  correction  before  the  discontinuity 
being  given  by  WWVTC  from  the  previous  card.  Note  that  if  the  discontinuity 
occurs  before  epoch,  i.e.  JMPDAY  MJDOCH,  then  the  correction  applied  at  the 
instant  of  epoch  will  be  WVArrC2. 

8,8  Observation  cards 

PROP  is  flexible  with  regard  to  the  variety  of  observation  formats  it  is 
capable  of  accepting.  It  can  deal  with  the  standard  formats  of  such  organiza¬ 
tions  as  the  Smithsonian  Astrophysical  Observatory  (SAO),  the  National  Aero¬ 
nautics  and  .Space  Administration  (NASA)  and  the  Radio  and  Space  Research 
Station  (RSRS).  It  can  also  be  very  easily  augmented  to  cope  with  new  formats, 
by  adding  extra  sections  to  the  basic  observations-input  subroutine  OBSIN. 
(Sections  can  equally  easily  be  removed,  for  computer  storage  economy,  if  they 
are  not  required.) 

Details  of  the  formats  covered  by  PR0P5  are  given  in  Appendix  B.  Here 
we  shall  merely  explain  how  a  complete  set  of  observation  cards  is  made  up. 

The  subroutine  OBSIN  reads  an  observations-format  or  ”LOOXEE**  card,  so 
called  from  the  name  of  the  corresponding  Fortran  variable.  This  card  causes 
control  to  be  directed  to  a  section  of  OBSIN  which  reads  observations  in  a 
specific  format.  The  several  formats  available  do  not,  in  general,  correspond 
to  particular  types  of  observation,  but  rather  to  the  styles  which  various 
organizations  have  adopted  in  publishing  satellite  observations.  So,  for 
example,  the  Moonwatch  format  for  presenting  observations  in  right  ascension 
and  declination  is  quite  different  from  the  RSRS  format. 

Each  such  section  of  OBSIN  reads  observation  cards  until  a  blank  card 
is  encountered,  after  which  OBSIN  interprets  the  next  card  as  a  new  LOOXEE 
card.  A  complete  set  of  observation  cards  for  a  given  run  will  thus  contain 
as  many  subsets  as  desired;  each  subset  consists  of  a  LOOXEE  card  followed  by 
observations  in  the  corresponding  format  and  concluded  with  a  blank  card. 

Exit  from  OBSIN  occurs  when  a  blank  LOOXEE  card  is  read.  This  means  that  the 
last  observation  of  the  complete  set  must  be  followed  by  two  blank  cards;  the 
first  of  these  concludes  the  last  subset,  and  the  second  concludes  the  set 
as  a  whole.  These  two  blank  cards  must  be  present  even  if  all  the  observations 
were  read  under  the  same  LOOXEE. 


Each  LOOXEE  card  contains  a  single  quantity,  read  under  format  (l6),  as 
follows : 

LOOXEE  (cols .1-6)  This  specifies  the  format  under  which  succeeding  cards  are 

to  be  read,  as  explained  above  and  detailed  in  Appendix  B. 

9  DESCRIPTION  OF  OUTPUT 

The  output  of  PROP  has  been  designed  to  be  as  nearly  self -exp  lematory 
as  possible.  The  purpose  of  this  section  is  to  make  such  comments  as  seem 
necessary  upon  the  output  of  a  typical  run*,  reproduced  as  Fig. 5.  This  is  the 
output  of  the  run  for  which  Fig. 4  lists  the  input  data. 

9.1  Preliminary  output 

The  first  three  computer  sheets  give  certain  preliminary  information 
relating  to  the  whole  run. 

The  title  is  followed,  in  brackets,  by  the  date  and  time  at  which  the 
computer  started  processing  the  job.  This  is  a  non-standard  facility,  added 
to  the  ICL  1907  version  of  PROP,  which  has  been  found  useful  in  the  day-to-day 
organization  of  orbit  determinations  on  the  R.A.E.  computer. 

No  comment  is  needed  on  the  listing,  for  reference,  of  the  control  para¬ 
meters  and  of  the  geophysical  constants  in  current  use.  Note  the  remark  which 
draws  attention  to  the  total  n\imber  of  zonal  harmonics  actually  to  be  used  (as 
specified  on  the  control  card).  In  fact  a  remark  appears  whenever  any  of  the 
constants  L,  EJ22  and  DENSCH  (S8.2)  is  reset. 

Printing  of  the  epoch  dace  and  satellite  identification  is  preceded, 
in  the  case  of  prediction  (S6.2),  by  a  listing  of  the  orbital  elements  used 
for  the  prediction. 

The  listing  of  sensor  locations  is  straightforward.  Note  that  the 
stations  numbered  1  to  18  (which  are,  in  fact,  NASA  stations)  have 
MARKER  =  2,  indicating  that  the  input  coordinates  were  X,  Y,  Z.  The  remain¬ 
ing  stations  have  MARKER  =  0,  indicating  the  the  input  coordinates  were 
latitude,  longitude,  height. 

The  listing  of  observations  needs  a  few  comments.  Each  observation  is 
assigned  a  serial  number  (in  the  example  they  run  from  1  to  57),  which  is  used 
for  reference  within  the  run  only.  If  an  observation  card  is  rejected  due  to 


*In  fact  this  run  is  not  entirely  "typical".  The  visual  observations  are 
of  poorer  quality  than  usual,  possibly  because  the  observers  were  concentrating 
on  flash  assessment. 
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an  anticipated  error  (i.e.  one  which  does  not  cause  deletion  of  the  Job), 
then  an  error  message  is  printed  by  the  progra,n,  the  information  from  that 
card  is  not  stored,  and  no  serial  number  is  assigned  to  it. 

The  type  number  refers  to  one  of  the  types  listed  in  section  2.2;  3 
indicates  direction-cosines  and  2  indicates  R, A. /declination  observations; 
if  R.A.  and  declination  have  been  converted  from  azimuth  and  elevation  this 
is  indicated  by  an  asterisk  before  the  type  number.  Note  that  information 
from  a  pair  of  direction-cosine  cards  (see  Fig. 4)  read  under  LOOXEE 
(see  Appendix  B)  is  printed  on  a  single  line. 

The  information  given  at  the  end  of  the  line  for  every  observation,  e.g. 
-27  REV  +  .27  (.41),  gives  a  rough  measure  of  the  spread  of  the  observations. 
The  quantity  -27  REV  +  .27  gives  the  difference  between  the  time  of  the  epoch 
for  the  run  and  the  time  of  the  paurticular  observation,  expressed  in  complete 
revolutions  plus  a  positive  fraction  of  a  revolution.  The  fraction  in  brackets, 
e.g.  (.41),  is  the  argument  of  latitude  (i.e.  argument  of  perigee  plus  true 
anomaly)  of  the  satellite's  position  at  the  time  of  the  observation. 

The  five  column  headings  from  OBS.  NO.  to  TIME  relate  to  all  the  observa¬ 
tions.  Other  column  headings  are  generated  each  time  a  LOOXEE  card  is  read, 
and  refer  to  a  particular  type  of  observation.  If  a  group  of  observations  do 
not  have  Individual  accuracies  assigned  to  them  but  are  given  standard  weights  - 
e.g.  with  direction-cosines  read  under  LOOXEE  444444  -  then  the  appropriate 
information  is  given  after  all  the  observations  in  that  group. 

The  transformation  of  direction-cosines  to  azimuth  and  elevation  is  done 
only  for  the  purpose  of  having  a  printed  record  of  the  observation  in  this 
form;  within  the  computer  the  observation  is  used  in  its  original  direction- 
cosines  form.  The  rule  that  obseived  quantities  are  always  used  in  the  form 
in  which  they  are  given  does  not  hold  for  azimuth/elevation  observations;  these 
are  converted  to  R.A. /declination,  and  the  computer  holds  a  record  of  the 
original  form  only  until  'processing  of  observations*  (S7)  is  complete. 

With  angular  observations  both  the  right  ascension  and  declination  and 
the  azimuth  and  elevation  are  listed.  Note  that  the  listing  is  given  by  sub¬ 
routine  OBSIN,  so  that  the  processing,  by  subroutine  PROCES,  has  not  yet  been 
incorporated,  and  the  conversion  of  angles  from  one  form  to  the  other  is  carried 
out  relative  to  the  standard  PROP  axes  (S5).  This  means  that  if  an  observation 
starts  as  right  ascension  and  declination  the  listed  azimuth  and  elevation  are 
only  approximate,  particularly  if  the  observation  relates  to  the  equinox  of 
1855,  1975  or  1900. 
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The  angular  accuracy  la  given,  In  seconds  of  sore,  under  S  (DEC).  This 
la  a  compounded  value  containing  a  component  derived  from  the  timing  error 
(S2.5).  Moonvatch  observations  read  under  LOOXEE  555555  do  not  have 
Individual  accuracies  attached  to  them,  and  they  are  all  assigned  a  standard 
a  priori  standard  deviation  by  the  program  (S2.5). 

The  preliminary  output  Is  concluded  by  a  listing  of  the  initial  estimates 
of  elements,  whether  these  were  obtained  directly  as  input  data  or  indirectly 
from  prediction.  The  elements  are  laid  ov.t  in  standard  tabular  form  correspoi  d- 
ing  to  the  •five-card  orbital  elements*  (S8.4.1),  with  the  addition  that  the 
derived  semi -major  aocis  (in  km)  is  printed  before  the  line  contai^.ing  the 
e-elements . 

N.B.  This  listing  of  initial  elements,  though  conveniently  thought  of 
as  concluding  the  preliminary  output,  is  actually  provided  at  the  beginning  of 
the  differential-correction  subroutine. 

9.2  Output  during  the  differential-correction  process 

Each  iteration  of  the  process  is  numbered.  If  the  mode  of  refinement 
(S6.4)  has  changed  prior  to  the  current  iteration,  a  statement  of  the  new  mode 
is  printed  at  the  head  of  the  table  of  residuals.  The  residuals  for  each 
observation,  both  the  weighted  and  the  unweighted  ones,  are  printed,  and  the 
units  of  the  unweighted  residuals  are  specified.  For  observations  of  right 
ascension  and  declination  a  further  quantity  is  printed;  this  is  a  "total" 

2  2  i 

unweighted  residual  given  by  {  (cos  6  fki)  +  (Al> )  where  ^  and  ^  are 

the  residuals  in  right  ascension  and  declination. 

The  marking  of  each  observation  with  three,  two,  one  or  no  asterisks 
relates  to  the  separating  of  observations  into  the  four  •rejection  levels* 
(s6.7).  Those  marked  with  three  asterisks  will  invariably  be  rejected  at  the 
end  of  the  iteration,  and  those  with  no  asterisks  will  be  accepted.  The 
acceptance  or  rejection  of  observations  marked  with  one  or  two  asterisks  depends 
on  rather  complicated  criteria  which  need  not  be  set  out  here;  the  program 
always  prints  a  comment  as  to  which  levels  are  being  rejected.  Note  that 
•residual*  has  the  sense  •observed*  minus  ‘computed* , 

The  text  at  the  end  of  each  iteration  is  self-explanatory .  Note,  however, 
that,  although  reference  is  made  to  the  number  of  "fresh  rejections",  there  is 
no  information  about  "fresh  acceptances",  i.e.  about  observations  rejected  on 
the  previous  iteration  but  now  accepted.  Prediction  of  the  value  of  e  after 
the  next  iteration  is  based  on  assumption  of  linearity. 
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9.3  Concluding  output 

Once  convergence  has  been  achieved,  standard  deviations  for  the  final  set 
of  elements  are  printed,  laid  out  in  the  same  way  as  the  corresponding  orbital 
elements.  If  there  are  orbital  elements  which  are  not  parar.eters  then  these 
have  not  varied  during  the  refinement  process  and  the  corresponding  numbers  in 
the  standard-deviations  layout  are  zero. 


The  symmetric  matrix  which  is  the  final  block  of  data  printed  consists  of 
the  correlations  of  the  orbital  parameters  in  the  order  (of  those  present ):- 


®2’ 


i 


n 


(j 


M 


At  this  final  stage  there  is  also  some  punching  of  cards.  Final  elements 

are  punched  on  five  cards ;  these  elements  differ  from  the  printed  set  in  being 

exclusive  (S4.3)  elements.  They  constitute  the  definitive  set  of  elements  for 

the  epoch.  The  covariance  matrix  of  the  orbital  parameters  is  punched  on 

(or  21  or  3'*  (‘'•'  .5))  cards.  Note  that  the  punched  covariance  matrix,  unlike 

the  printed  correlations,  refers  to  non-standard  parameters  if  these  have  been 

used  during  the  differential  correction  (S6.5).  Also  note  that  the  units  are 

those  used  inside  the  computer  (S5);  e.g.  the  printed  standard  deviation  for 

-2C 

M.  is  0.00057®/ day,  but  the  punched-card  variance  is  1.55  X  lO 
2 

(radians/sec)  . 


Output  in  the  event  of  failure 

Program  failures  may  be  considered  to  be  of  three  types.  (Rejection  of 
an  observation  card  due  to  an  anticipated  error  (S9.1)  is  not  counted  as  a 
fai lure. ) 


First  there  are  the  standard  failures  referred  to  in  section  6.1  .  A 
failure  of  this  type  can  only  occur  after  input  of  all  the  data  for  a  given 
run.  In  PR0P3  there  are  six  possible  standard  failures  (the  first  occurs  in 
the  main  program  and  the  others  in  subroutine  DIFCOR)  and  the  corresponding 
output  is  as  follows 


DISCONTINUE  -  TOO  FEW  OBSERVATIONS  INPUT  AND  ACCEPTED 
DISCONTINUE  -  TOO  FEW  OBSERVATIONS  LEFT 

DISCONTINUE  -  ELEMENT  TYPE  INCONSISTENT  WITH  USE  OF  PREVIOUSLY- DETERMINED  ORBIT 

DISCONTINUE  -  DEGREES  OF  FREEDOM  NEGATIVE 

DISCONTINUE  -  LIMIT  ON  ALLOWED  NUMBER  OF  ITERATIONS  REACHED 

DISCONTINUE  -  DIVERGENCE  ON  TWO  SUCCESSIVE  ITERATIONS. 
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Second  there  are  failures  which,  though  they  lead  to  deletion  of  the 
program  (via  a  Fortran  STOP  instruction),  have  been  anticipated  and  provide 
explanatory  printing  before  the  deletion.  These  all  occur  (or  may  occur) 
before  input  of  the  complete  data  for  the  run,  and  this  is  why  it  has  not 
been  arranged  for  the  next  run  to  be  started. 

Finally  there  are  unanticipated  failures  which  will  leave  no  direct 
explanation  of  their  cause.  These  will  arise  from  breaking  a  fundamental 
rule  of  the  computer,  e.g.  through  overflow  or  a  data  card  punched 
inconsistently  with  the  associated  format  statement,  and  lead  to  deletion  of 
the  program  via  the  computer *s  supervisor  or  executive. 

10  PRESENTATION  OF  RESULTS 

It  will  be  normal  for  the  PROPJ  user  to  present  his  accumulated  results 
in  tabular  form,  particularly  wnen  a  number  of  orbit  determinations  for  a 
rivpn  satellite  have  been  carried  out  at  regular  inter'’’als.  It  may  also  be 
appropriate  to  illustrate  the  variation  of  certain  orbital  elements,  for 
example  the  eccentricity,  in  graphical  form,  but  this  section  considers  only 
the  preparation  of  tables. 

All  PR0P3  orbit  determinations  relate  to  midnight  epochs,  and  an  obvious 
first  quantity  to  tabulate  is  the  date  of  epoch.  It  is  a  good  practice  to  give 
this  both  as  a  noimal  date  and  as  an  MJD  number.  The  remaining  quamtities  to 
be  tabulated  will  normally  include  the  orbital  parameters,  with  their  standard 
deviations,  some  useful  derived  parameters  (functions  of  the  official  para¬ 
meters)  such  as  semi -major  axis  and  perigee  height,  also  with  standard 
deviations,  and  other  quantities,  such  as  the  number  of  observations  used. 

It  is  suggested  that  values  of  the  orbital  parameters  be  rounded  such 
that  standard  deviations,  when  similarly  rounded,  are  given  to  ore  significant 
figure,  or  at  most  two.  Thus  it  is  absurd  to  refer  to  an  eccentricity  value 
of  0.018^25  ±0.001264;  with  an  obvious  space-saving  convention  ore  would  round 
to  0.0184  IJ.  However  this  principle  raises  problems  when  certai:.  of  the 
parameters  are  highly  correlated.  There  are  two  possible  reasons  for  high 
correlations:  one  or  more  of  the  parameters  are  inherentlj^  il] -defined  due  to 
proximity  to  one  of  the  singularities  in  e  and  i  (s6.5);  or  else  the  cover¬ 
age  of  the  orbit  by  the  available  observations  is  poor. 

To  exemplify/  the  problem  just  referred  to,  consider  a  satellite  for 
which  e^  -  '.ool  ,  A  given  orbit  determination  may  yield  standard  deviations 


20C, 


k7 


including  o(fi^)  -  0®.005  and  o(w^)  =  o(m^)  =  2*;  i.e.  is  worth  quoting 

to  three  decimal  places,  but  w  and  M  onlj'  to  the  nearest  degree.  However, 

o  0 

and  will  have  a  correlation  verj'  close  to  -1  .0,  and  the  composite 

parameter  M  w  may  be  loiown  as  accurately  as  fJ  .  To  conform  to  the 

0  0  o 

general  principle  of  the  last  paragraph,  while  retaining  real  accuracy,  the 
PROP  user  should  quote  to  the  nearest  degree  and  to  three 

decimal  places,  not  quoting  itself  at  all  (-  there  is  no  problem  here  with 

Mj  etc.).  The  standard  deviation  of  the  composite  parameter  would  have  to  be 
quoted  instead  of  o(M^),  of  course,  and  this  is  given  by 

2  2  2 
o(m+w)  =o(a))+r  coy(w  ,  M  )  +  0  (m  ) 

'o  o  'o  'o  o  'o 


An  attempt  tc  calculate  +  “q)  using  the  value  of  cov(w^»  M^)  on  an 

output  covariance-matrix  card  will  fail,  since  these  cards  do  not  give  enough 

significant  figures  to  cover  the  high-correlation  situations,  but  this  should 

be  unnecessary  anyway:  the  user  should  have  chosen  a  JELTYP  value  (sS.2)  such 
2 

that  o"(Mq  is  given  directly  (s6.5). 

It  is  suggested  that  M  w  be  quoted,  instead  of  M  >  whenever  a 

0  0  0 

satellite’s  eccentricity  is  less  than  about  0.05.  This  is  only  a  guide,  of 

course;  one  should  not  change  parameters  from  line  to  line  of  a  table,  just 

because  e  is  oscillating  about  0.05. 

When  i  is  near  0®  or  1^0®  (but  assuming,  first,  that  e  is  not 
0  o 

small)  there  is  a  similar  problem  for  (M^  being  now  all  right)  and 

u)  +  fi  or  w  -  n  should  be  quoted  instead  of  <*>  .  It  will  be  recalled 
O  O  O  0  o 

(S6.5)  that,  with  the  appropriate  choice  of  JELTYP,  the  covariance  cards  will 

2 

give  0  (u)  -t  k  n  ),  but  this  will  certainly  give  o(w  -t  fi  )  or  o(i'  -  f!  ) 

0  0  o  0  o  o 

to  sufficient  accuracy,  since  k  will  be  near  to  +1  or  -1.  When,  in  addition, 

e  is  small  it  will  be  necessary  to  quote  parameters  fi  ,  u)  4  fi  smd 
o  0  0  o 

M  +  w  +  fi  ,  or  else  fi  ,  u  -  fi  and  M  -*■  <*)  -  ;  no  new  principle  is 

000  000  OOQ  r 

involved . 

Returning  to  general  aivi.  e  on  the  quantities  to  tabulate,  e^,  i^  etc. 
will  of  course  be  tabulated  if  they  are  parameters.  If  they  are  merely  elements, 
for  which  the  exclusive  values  (S4.5)  are  zero,  it  may  be  preferred  to  omit 
them.  If  e^,  i^  etc.  are  quoted  it  will  certainly  be  necessary  to  remind 
readers  of  the  difference  between  inclusive  and  exclusive  values. 
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Other  quantities  which  it  may  be  desired  to  tabulate  are:-  number  of 
observations  d,  counting  any  rejected  during  ainalysis;  number  of  observa¬ 
tions  rejected;  n\unber  of  days  spanned  by  the  observations  (e.g.  ^.8);  and 
the  final  value  of  e  (e.g.  1.8). 

The  best  way  to  prepare  am  error-free  table  is  to  use  the  output  cards 
from  each  run.  A  Fortran  program  TOP  (-  Tabulation  of  Orbital  Parameters) 
has  been  written  which  reads  in  as  many  sets  of  these  as  desired  and  prints 
the  table,  one  row  for  each  set  of  output  cards.  As  the  output  cards  cover 
only  elements  and  covariances,  an  extra  card  is  required  with  each  set,  giving 
the  number  of  observations,  etc.  TOP  has  been  written  with  formats  suitable 
for  the  listing  of  orbital  parameters  for  Ariel  5,  as  determined  from  Mini¬ 
track  observations,  but  it  can  easily  be  modified  for  use  with  other 
satellites . 

11  CONCLUSIONS 

We  have  described  in  detail  the  operation  of  PROP,  a  Fortran  program 
for  refining  the  orbital  parameters  of  an  earth-satellite.  Although  the 
program  is  not  capable  of  refining  station  coordinates  or  geophysical  constants, 
its  flexibility  in  other  ways  makes  it  an  excellent  tool  in  the  analysis  of 
the  orbits  of  many  close  earth-satellites.  However,  certain  limitations  in  the 
dynamical  model  should  be  kept  in  mind,  to  avoid  using  PROP  for  tasks  for 
which  it  is  unsuited. 

(i)  The  program  does  not  contain  formulae  for  representing  luni-solar 
perturbations.  This  fact  imposes  an  upper  limit  on  the  value  of  a  satellite's 
semi-major  axis,  beyond  which  the  program  should  not  be  used.  (As  a  very 
rough  rule-of-thumb  this  upper  limit  might  be  taken  to  be  two  earth-radii.) 

Some  allowance  con  be  made  for  luni-solar  effects  by  calculating  the  main 
secular  terms  in  e,  i,  fi  and  w  before  using  PROP;  these  terms  cem  then  be 
fed  in  as  fixed  elements  in  the  model.  It  is  planned  to  include  this  com¬ 
putation  internally  in  the  next  version  of  PROP. 

(ii)  The  representation  of  atmospheric-drag  effects  is  not  adequate 
for  a  satellite  strongly  affected  by  drag.  This  imposes  a  lower  limit  on  the 
value  of  a  satellite's  perigee  distance,  below  whi  h  PROP  should  noi  be  used. 
This  lower  limit  depends  on  so  many  factors  -  li  •  iding  the  satellite’r  drag 
coefficient  and  the  level  of  solar  activity  -  that  no  clear  cri  ter  ion  an  be 
suggested  for  assessing  it;  the  best  guide  ir  ; robaP  thp  value  of  M^. 
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In  the  next  version  of  PROP  it  is  planned  to  deal  with  the  hlgh-drag 
situation  by  a  standard  procediire:  to  alter  the  weighting  of  the  obseinrations 
used  in  a  given  orbit  determination  so  that  the  along-track  'component  of  each 
observation  (the  component  most  affected  by  drag)  receives  less  weight  than 
the  cross-track  components . 

(iii)  The  effects  of  all  tesseral  harmonics  have  been  neglected,  with 
the  exception  of  the  dominant  J.  *  effect.  This  leads  to  errors  of  the  order 
of  ^  km  in  the  position  of  a  satellite  as  computed  from  a  given  set  of  orbital 
elements.  It  is  hoped  to  improve  this  situation  in  the  next  version  of  PROP. 

The  authors  have  attempted  to  be  comprehensive  in  preparing  this  manual, 
but  it  is  unlikely  that  there  are  no  errors  of  omission  or  commission  in  the 
text.  In  anticipation  of  the  probability  that  it  will  be  revised  and  reissued 
when  PR0P4  is  available  we  invite  all  readers  of  the  document,  and  in  particular 
those  who  actually  use  the  program,  to  bring  any  queries  or  criticisms  to  our 
attention. 
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BLOCK  DATA 
Main  prog. 

^^\cos 

ALTELS 

AZEL 

CONVER 

CONVRT 

COTOEL 

DHMS 

DIFCOR 

EAFKEP 

(b ) 

^  '^ECENTR 
ELPRIN 
ELREAD 
E PARAS 
GEOCOR 

WolA 

MODJD 

riEVfELS 

NUPOLY 

OBSIN 

PARSHL 

PRDATE 

PRELON 

PRENUT 

PROCES 

PVREAD 

REFCOR 

REJECT 

SATELS 

SATXYZ 

SDPRIN 

sense; 

SHOPER 


Sets  certain  constants. 

Sets  other  constants  and  supervises  flow  -  MASTER  PROP  in  1907 
Fortran. 

Arc-cosine . 

Alters  orbital  elements  (‘'exclusive”  to  and  from  "inclusive"). 
Az/el/range  from  topocentric  satellite  position. 

Tests  refinement  process  for  convergence. 

Converts  between  az/el,  R.A./dec  and  direction  cosines. 
Coordinates  to  osculating  elements. 

Given  days,  or  degrees,  introduces  hours,  minutes  and  seconds. 
Differential  correction  -  the  refinement  process. 

Eccentric  anomaly  from  Kepler's  equation. 

Geocentric-to-topocentric  conversion  of  satellite  position  etc. 
Prints  (and/or  punches)  orbital  elements. 

Reads  orbital  elements. 

Abstracts  list  of  parameters  from  the  array  of  elements. 
Geographical  coordinates  from  geocentric  Cartesians. 

Matrix  division  or  inversion  (single-length). 

Expresses  date  as  MJD  number. 

Sets  new  values  of  parameters  after  each  refinement. 

Nearest  integer. 

New  polynomial  coefficients  after  shift  of  origin. 

Input  of  observations  in  various  formats. 

Partial  derivatives  of  (geocentric)  satellite  coordinates. 

Prints  date,  given  MID  number. 

Sets  time-independent  quantities  for  perturbations. 

Correction  of  R.A./dec  observations  for  precession  and  nutation. 
Processing  of  observations  prior  to  refinement  processes. 

Reads  position  and  velocity. 

Refraction  correction  of  observed  elevation. 

Decision  as  to  observations  to  be  rejected. 

Position  and  velocity  from  mean  orbital  elements. 

Perturbed  geocentric  position  and  velocity. 

Prints  standard  deviations  of  elements. 

Input  of  sensor  locations. 

Computes  short-periodic  perturbations. 
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SIDANG 

SINXOX 

TRINV 

XLONG 


Sidereal  angle  relative  to  1950.0. 

Sin  X  /  X,  (1  -  cos  x)  /  X  and  (x  -  sin  x)  /  x^. 
Polar  coordinates  from  Cartesians  (two-dimensional). 
Computes  long-periodic  perturbations. 


Notes 

(a)  The  ACOS  and  NINT  functions,  though  not  included  in  A.S.A.  Fortran, 
are  provided  by  both  the  Atlas  and  the  1907  computers  and  need  not  be 
provided  separately  for  these  computers. 

(b)  ECENTR  is  called  as  the  dumn^-argument  subroutine  CENTRE,  this  being 
organized  by  an  EXTERNAL  statement  in  the  main  program.  For  orbiting  bodies 
other  than  earth  satellites  a  different  CENTRE  may  be  used. 

(c)  MAOlA,  similarly,  is  called  as  MAOI  .  The  *A*  version  gives  single¬ 
length  matrix  inversion.  If  double-length  inversion  is  required,  MAOlB 
should  be  used  instead,  and  references  to  MAOlA  in  the  main  program  replaced 
by  references  to  MAOlB, 
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Appendix  B 

OB3ERVATIONS-FORMATS  CURRENTLY  PROVIDED  FOR  BY  THE  SUBROUTINE  OBSIN 

Each  observations-format  is  associated  with  sin  appropriate  identifier. 
Before  any  group  of  observations  is  read,  the  appropriate  identifier  is 
read  into  the  Fortran  variable  LOOXEIE,  and  it  is  conv  enient  to  list  the 
formats  currently  available  under  the  corresponding  LOOXEE  numbers, 

B.l  LOOXEE  mm 

0 

This  format  reads  Baker-Nunn  observations  in  the  standard  SAO  format  , 
implemented*  in  OBSIN  by : 

FORMAT  (17,  5X,  15,  I6,  212,  F6.i*,  1X,  212,  F5.5,  Al,  212,  F4.2,  II,  12,  IX, 
n,  7X,  T7.k)  . 

The  17  quantities  input  are  read  into  the  following  Fortran  variables: 

This  is  the  satellite  identification  according  to  the 
international  designation;  e.g.  satellite  1967“^9b  is 
punched  as  6704902. 

This  is  the  number  of  the  observing-station;  e.g.  Maui, 
Hawaii,  is  b90l2. 

This  is  the  date  in  6-digit  form;  e.g.  1963  May  22 
is  680522. 

These  give  the  time  of  the  observation;  e.g.  IJh  4m 
29.01-''>  is  punched  as  150429651b.  (Note  that  column 
55  must  be  blank  unless  the  time  is  given  to  0.1  msec.) 
These  give  the  satellite's  right  ascension,  e.g. 

5h  44m  19s  is  punched  as  b544l9bbb.  A  decimal  point 
is  implicit  between  columns  40  and  41  . 

These  give  the  satellite's  declination,  a  decimal 
point  being  implicit  between  columns  50  and  51;  e.g. 

_ro  291  52". 41  is  punched  as  -05295241  and 
*5°  29'  52”. 41  is  punched  as  bb5295241  .  The  minus  sign 
(if  any)  must  appear  in  column  44. 


JDENT  (cols  .1-7) 

NOSTAT  (cols .15-1 7) 

MDAY  (cols. 19-25) 

MliR  (col 3. 24-25 
MIN  (cols. 26-27)  [ 

OEC  (cols. 29-^5)  I 
>-WAH  (cols. 55-56) 
MRAM  ( cols. 57-59 ) ' 
RAE  (col  3. 59-45) 

MIGI.I  (-01.44) 
MDECD  (cols. 45-4^ 
MOF.^M  (''ol  3 .47-49) 
n'.CG  (cols. 49-52) 


“Only  S.A.O.  R.A./dec  observations  are  covered  at  jiresent.  It  is  hoped 
t  1  f'ovr>r  other  C.A.O.  formats  {i/m,  A/E  etc.)  in  an  extended  LOOXEE  111111  in 
n  later  version  of  PROP. 
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NTSD  (col. 55)  This  is  a  code  digit  indicating  the  timing  accuracy  of 

g 

the  observation  ;  zero  (or  blank)  signifies  that  there 
is  no  estimate  of  timing  accuracy,  in  which  case  the 

g 

subroutine  sets  an  arbitrary  accuracy  of  0  .005;  T 

g 

signifies  an  accuracy  of  0  .0001 ;  2  signifies  an  accuracy 
of  O^.OOl;  5  signifies  an  accuracy  of  0^.005;  and  4 

g 

signifies  an  accuracy  of  0  .01  *. 

NASD  (cols. 54-55)  This  is  the  angular  accuracy  of  the  observation  in  seconds 

of  arc. 

INECiEQ  (col. 57)  This  is  a  code  digit  indicating  the  date  of  the  equator 

and  equinox  to  which  the  observations  are  referred;  the 
digits  0,  1,  2,  3  and  4  refer  respectively  to  the  equator 
and  equinox  of  date,  1‘^j55.0,  1^75.0,  1 900.0  Hpd  1950.0. 
TIMOR  (cols. 65-71)  This  is  the  time  correction  to  be  subtracted  from  the  time 

in  the  Al  system  to  give  time  in  UTl  ;  e.g.  a  difference 
of  42"’. 21 7  is  punched  as  b422l7b.  The  decimal  point  is 
implicit  between  colvimns  67  and  60. 

The  card  may  contain  other  quantities  ,  e.g.  a  reference  number  in  cols, 8-1 2, 
which  are  not  used  by  PROP. 

B.2  LOOXEE  222222 

This  format  reads  direction-cosine  (Minitrack)  or  reinge  or  range-rate 
observations  in  a  particular  form.  At  R.A.E.  it  has  been  used  only  for  the 
satellite  OGO-2  and  may  be  abandoned  in  future  versions  of  PROP.  The  format 
is  implemented  in  OBSIN  by: 

FORMAT  (12.  0X,  15,  F9.9,  IOX,  II ,  2(4x,  II ,  El  5.0)) 

The  IO  quantities  input  are  read  into  the  following  Fortran  variables: 
NOSTAT  (cols .1-2)  This  is  the  number  of  the  observing-station .  Note  that 

this  is  a  2-digit  code,  covering  a  possible  maximum  of 
99  stations,  and  is  inconsistent  with  the  international 
4-digit  code  for  station  numbers.  (When  used  for  OGO-2 
the  observations  were  from  the  following  stations, 
numbered  (j1  to  1 5 :  College,  Blossom  Point,  Fort  Myers, 

Grand  Forks,  Johannesburg,  Lima,  Mojave,  Newfoundland, 
Woomera,  Quito,  Santiago,  Winkfield,  Rosman  (u.S.A.), 
Tananarive  (Malagasy  Republic)  and  Carnarvon  (Australia).) 

*Ref.0  also  covers  values  of  4  to  0,  but  these  are  not  :overed  in  PR0P3. 
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NTIMHR  (cols. 11-15)  The  number  of  complete  hours  elapsed  since  MJD  56955.0 

(i960  January  0.0)  to  the  time  of  the  obseirvatlon;  e.g. 
1200  hours  on  I965  November  9  (viz.  MJD  59075  +  12  hours) 
is  2140  days  and  12  hours  from  MJD  56955.0,  so  it  is 
punched  as  51572. 

TIM  (cols .16-24)  The  fraction  of  an  hour  elapsed  from  NTIMHR  to  the  tlmf 

of  the  obseivation;  e.g.  27m  45.27s  is  punched  as 
46201 9444.  The  decimal  point  is  implicit  before  coluinn 

16. 

NOBS  (col. 55)  The  number  of  observed  quantities  given  on  the  card:  1 

for  a  range  observation  or  a  range-rate  observation,  2 
for  a  direction-cosines  observation. 

NTYPl  (col.40)  A  code  digit  indicating  the  type  of  observed  quantity 

which  is  read  into  OBSl  :  5,  4  and  7  refer  respectively 
to  range,  a  range-rate  and  an  east-west  direction  cosine. 

OBSl  (cols. 41 -55)  The  first  observed  quantity;  rsmge  must  be  in  units  of 

half -earth-radii,  where  the  earth's  radius  is  taken  to  be 
6378.165  km;  range-rate  must  be  in  units  of  earth-radii/ 
hour;  direction-cosines  are,  of  course,  dimensionless. 

A  decimal  point  and/or  exponent  should  be  present. 

NTYP2  (col. 60)  A  code  digit  Indicating  the  type  of  quantity  read  into 

OBS2;  only  one  value  is  recognised:  8  for  a  north-south 
direction  cosine. 

0BS2  (cols. 61 -75)  The  north-south  direction  cosine. 

Note:  If  NOBS  =  1,  indicating  a  range  or  range-rate  obseivation,  the  fields  of 
NTYP2  and  0BS2  will  be  blank. 

B.5  LOOXEE  555555 

This  format  reads  range  or  range-rate  observations  in  the  standard  form 

supplied  by  NASA.  It  is  implemented  ^n  OdSIN  by: 

FORMAT  (7X,  2A5,  IX,  I6,  IX,  212,  IX,  f6.4,  loX,  F9.5,  F7.5) 

The  7  quantities  input  are  read  into  the  following  Fortran  variables: 

NAMEST(2)  (cols. 8-1 5)  The  station  name,  in  a  6-digit  alphanumeric  code.  NAMEST 

must  Identify  with  one  of  the  station  names  held  (as 
text)  in  the  array  M0NICA(56).  Each  of  these  names  is 
held  in  two  words  of  storage;  thus  the  first  name  is  held 
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in  M0NICA(1  )  and  M0NICA(2),  and  the  name  is  held  in 

th 

M0NICA(2n-1)  and  M0NICA(2n).  If  NAMEST  is  the  n 
name,  station  coordinates  will  be  taken  from  a  sensor 
card  for  which  the  station  number  is  n. 

The  18  code-names  currently  held  in  MONICA  are, 
in  order,  COLEGf,  BP0IN6,  FTMYR6,  GF0RK6,  J0BUR6, 

LIMAP6,  M0JAV6.  NEWFL6,  00MER6,  QUITO6,  smG6, 

WKKFL6,  ROSRAN.  mADGAR,  CARVON,  ULASK6,  MADGA6  and 
0R0RA6.  Thus,  for  example,  the  station  name  QUITOo 
refers  to  a  sensor-card  numbered  lO. 

MDAY  (cols .15-20)  The  date  in  6-digit  form;  e.g.  I968  May  22  is  punched  as 

680522. 

MHR  (cols .22-25) 'i  The  hours  and  minutes  components  of  the  UTC  time  of  the 

MIN  (cols .24-25)  J  observation;  e.g.  07h  28m  is  punched  as  b728. 

SEC  (cols. 27-52)  The  seconds  component  of  the  time  of  the  observation; 

2 

e.g.  17  .819  is  punched  as  17819b.  The  decimal  point 
is  implicit  between  columns  28  and  29. 

RANGE  (cols. 49-57)  The  satellite *s  range  from  the  station,  in  kilometres, 

the  decimal  point  implicit  between  columns  5^  and  55; 
e.g.  1936.565  km  is  pijuiched  as  bbl 956565. 

RGEDOT  (cols. 58-64)  The  rate  of  change  of  range,  in  km/sec,  the  decimal 

point  implicit  between  columns  59  and  60;  e.g.  4.14812 
km/sec  is  punched  as  b4148l2. 

I 

N.B.  Only  one  of  RANGE  and  RGEDOT  may  be  present  on  the  same  card.  The  other 
field  must  be  left  blank.  (PROP  reads  both  fields,  but  interprets  the 
observation  as  of  "range”  or  "range  rate"  according  to  whether  RANGE 
is  not,  or  is,  zero.) 

B.4  LOOXEE  444444 

This  format  reads  standard  Minitrack  direction-cosine  cards,  and  is 
implemented  in  OBSIN  by: 

FORMAT  (7X,  2A5,  1 X,  I6,  1 X,  212,  IX,  F5.5,  28x,  f8.6,  4x,  II  ) 

Pairs  of  cards  are  looked  for,  and  single  cards  are  rejected.  The 
quantities  on  the  first  card  of  a  pair  are  read  into  the  following  Fortran 
variables : 
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NAME?r(2)  (cols. 8-1 5)  The  station  name,  in  a  6-digit  alphanumeric  code,  as 

described  under  LOOXEE  35553?. 

MDAY  (cols. 15-20)  The  date,  in  6-digit  form;  e.g.  1968  May  7  is  punched 

as  680507. 

MHP  (cols. 22-25)  I  The  hours  and  minutes  components  of  the  UTC  time  of  the 

MUf  (cols. 24-25)  J  observation;  e.g.  I6h  55^  is  punched  as  I655. 

SEC  (cols.27-?l)  The  seconds  component  of  the  time  of  the  observation;  e.g. 

54  .43  is  punched  as  5^^5b.  The  decimal  point  is  implicit 
between  columns  28  emd  29. 

DCL  (cols. 60-67)  The  east-west  direction-cosine,  with  the  decimal  point 

implicit  between  columns  61  and  62. 

NP  (col. 72)  This  is  for  checking  that  DCL  is,  in  fact,  the  east- 

west  direction-cosine.  If  NP  ♦  2  the  check  fails  and  the 
card  is  rejected. 

The  first  four  quantities  on  the  second  card  of  the  pair  are  read  into 
the  Fortran  variables:  NAMEST,  MMY,  MHR,  MIN  exactly  as  before;  in  fact  these 
quantities  overwrite  those  which  were  set  from  the  first  card.  The  remaining 
quantities  are  read  into  the  following  Fortran  variables  :- 

SEC2  (cols  .27-51)  This  must  be  within  5  msec  of  SEC,  otherwise  the  pair  of 

data-cards  is  rejected. 

DCM  (cols. 60-67)  The  north-south  direction-cosine,  the  decimal  point 

Implicit  between  columns  61  and  62. 

NP  (col. 72)  This  is  for  checking  that  DCM  is,  in  fact,  the  north- 

south  direction-cosine.  If  NP  *  5.  the  check  fails; 
but  a  further  check  is  made:  if  NP  =  2  the  subroutine 
assumes  that  the  second  card  of  a  pair  has  been  lost. 

It  therefore  sets  DCL  =  DCM  and  SEC  -  SEC2  and  assumes 
that  this  card  is  the  first  of  a  pair.  If  NP  is 
neither  5  nor  2  the  card  is  rejected  (and  so,  of  course, 
is  the  previous  one  of  the  pair). 

Note:  NASA  Minitrack  cards  contain  a  5-digit  satellite  identification  (cols. 2- 
6),  but  this  is  not  used  by  PROP. 

B.5  LOOXEE  535555 

This  format  reads  angular  observations  on  standard  Moonwatch  cards  as 
supplied  by  the  S.A.O.  The  description  will  be  abbreviated  here,  since  the 
PROP  user  will  not  normally  be  punching  cards  in  this  format.  If  he  wishes  to 
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punch  optical  or  radar  observations,  the  format  described  in  section  B.6  below 
will  be  found  to  be  more  flexible  and  comprehensive,  and  its  use  is 
reconnnended. 

The  format  is  implemented  in  OBSIN  by: 

FORMAT  (17,  5X,  15,  I6,  212,  F6.i*,  15,  12,  F2.0,  5X,  Al,  212,  F2.0,  6X,  tl  )  , 


where  the  quantities  are  read  into  the  following  Fortran  variables: 


JDENT  (cols. 1-7) 
NOSTAT  (cols .13-1 7) 
MDAY  (cols  .18-23) 


MHR 

MIN 

SEC 

Ml 

M2 

AM3 


(cols. 24-53) 


(cols .34-40) 


Satellite's  international  designation. 

Observing  tstatlon's  number. 

Date  of  the  observation,  in  6-digit  form. 

Time  of  the  observation,  the  decimal  point  implicit 

g 

between  columns  29  eind  30;  e.g.  I9h  49m  33  .6  is  punched 
as  I949336bbb. 

Hours,  minutes  and  seconds  components  of  right  ascension, 
or  degrees,  minutes  and  seconds  components  of  azimuth. 
E.g.  an  R.A.  of  21h  50m  I2s  is  punched  as  b2l50l2,  an 
azimuth  of  115“  ^6*  18”  as  1154618. 


MSIGNl  ] 

^  [  (cols. 44-50) 

N2 

;ij43 

INEQEC,  (col. 57) 


Sign,  degrees,  minutes  and  seconds  components  of 
declination  or  of  elevation.  E.g.  a  declination  of 
-74®  17'  is  punched  as  -74l7bb,  an  elevation  of  44®  6*  14" 
as  b44o6l4. 

Code  digit  indicating  the  type  of  observation  and  the 
date  of  the  equator  eind  equinox  to  which  it  is  referred; 
zero  (or  blank)  implies  azimuth  and  elevation;  1,  2,  5 
and  4  imply  right  ascension  and  declination,  and  refer 
respectively  to  the  equator  and  equinox  of  1855.0,  1875.0, 
1900.0  and  1950.0. 


B.6  LOOXEE  666666 


This  format  reads  range  and/or  angular  observations  on  cards  pimched 
according  to  the  format  agreed  between  R.A.E.  (Famborough),  RSRS  (Slough), 
the  Royal  Greenwich  Observatory  ( Herstmonceux),  the  Royal  Observatory, 
Edinburgh,  and  the  Sub-Committee  on  Optical  Tracking  of  the  British  National 
Committee  for  Space  Research.  It  is  therefore  the  standard  British  format  for 
reporting  satellite  observations.  It  is  implemented  in  OBSIN  by: 
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FORMAT  (17,  512,  F6.i»,  F5.4,  1x,  II,  2F8.0,  f4.0,  I1,  f8.5,  F5.5)  , 


where  the  quantities  are  read  into  the  following  Fortran  variables : 


JDENT  (cols. 1-7)  The  satellite  designation  in  the  7-dlglt  international 

code;  e.g.  satellite  1967-^8b  is  punched  as  6704802, 
NOSTTAT  (cols, 8-11)  The  normal  4-diglt  number  assigned  to  a  station  by  a 

coordinating  centre  on  behalf  of  COSPAR;  e.g.  Hbrsebrldge 
is  5657. 


MY 

MM 

MD 

mhr' 

MIN 

SEC 


(cols  .12-17) 


(cols  .18-27) 


TSD  (cols. 28-32) 


NCODE  (col, 34) 


RAAZ  (cols. 35-42)  '1 

DECEL  (cols. 43-50)  / 


The  date  of  the  observation,  in  6-digit  form;  e.g.  1968 
May  7  is  680507. 

The  time  of  the  observation  (OTl);  e.g.  I6h  4m  6.254s 
is  punched  as  l6o4o6234b.  Note  that  the  decimal  point 
is  implicit  between  columns  25  and  24. 

The  timing  accuracy  in  seconds  of  time;  the  decimal 
point  is  Implicit  between  columns  28  and  29.  E.g,  0®.5 
is  punched  as  b5bbb. 

A  code  digit  which  indicates  the  format  and  units  of  the 
angular  part  of  the  observation  (see  below). 

The  first  and  second  position  coordinates,  l.e.  the  angular 
part  of  the  observation,  according  to  the  following  table: 


NCODE  value 

First  position  coordinate 

Second  position  coordinate 

0 

Right  ascension  measured  in 
hours,  minutes  and  seconds, 
a  decimal  point  implicit 
between  columns  40  and  41, 
e.g.  4h  I9ffl  58.4s  is  punched 
as  b41 9384b 

Declination  measured  in  degrees, 
minutes  and  seconds  of  arc,  a 
decimal  point  Implicit  between 
columns  49  and  50,  e.g.  -8* 

12*  7". 7  is  punched  as 

-0812077 

The  observation  has  been  corrected  for  both  components  (astro¬ 
nomical  and  parallactic)  of  refraction. 

1 

Format  Identical  to  that  for  NCODE  =  0,  but  the  observation  has 
been  corrected  for  astronomical  refraction  only. 

2 

Right  ascension  measured  in 
hours  and  minutes,  a  decimal 
point  implicit  between  columns 
38  and  39,  e.g.  4h  19.64m  is 
punched  .1^41 964^^ 

Declination  measured  in  degrees 
and  minutes  of  arc,  a  decimal 
point  iit5>licit  between  columns 

47  and  48,  e.g.  -8*  12*. 128  is 
punched  as_Q3,2128 

The  observation  is  corrected  tor  astronomical  refraction  only. 
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NCODE  value 

First  position  coordinate 

Second  position  coordinate 

5 

Right  ascension  measured  in 
hours  and  minutes,  as  with 
NCODE  -  2,  e.g.  4h  l9.f-^m  is 
punched  as 

b41 964bb 

Declination  measured  in 
degrees,  the  decimal  point 
implicit  between  columns  45 
and  46,  e.g.  -8®. 20233  is 
punched  as 

-0820253 

The  observation  has  been  corrected  for  astronomical  refraction 
only. 

h 

Azimuth  measured  in  degrees, 
minutes  and  seconds,  a 
decimal  point  implicit  between 
columns  41  and  42,  e.g. 

194*  5*  17". 4  is  punched  as 

19403174 

Elevation  measured  in  degrees, 
minutes  and  seconds,  a 
decimal  point  Implicit  between 
columns  49  and  50,  e.g. 

46®  19*  46". 6  is  punched  as 
b46l  9466 

The  observation  has  been  corrected  for  both  components  of 
refraction. 

5 

Azimuth  measured  in  degrees 
and  minutes,  a  decimal  point 
implicit  between  columns  39 
and  4o,  e.g.  194*  5*. 29  is 
punched  as 

1 940529b 

Elevation  measured  in  degrees 
and  minutes,  a  decimal  point 
implicit  between  columns  47 
and  48,  e.g.  46®  19.777  is 
punched  as 

b46l9777 

The  observation  has  been  corrected  for  both  components  of 
refraction. 

6 

Azimuth  measured  in  degrees, 
a  decimal  point  implicit 
between  columns  37  and  30, 
e.g.  194*. 05483  is  punched  as 

1 9405483 

Elevation  measured  in  degrees, 
a  decimal  point  implicit 
between  columns  45  and  46, 
e.g.  46®. 52965  is  punched  as 
b4652963 

The  observation  has  been  corrected  for  both  components  of 
refraction. 

7, 0,9 

Format  and  units  the  same  as  for  4,  5,  6,  respectively,  but 
the  observations  are  completely  uncorrected  for  refraction. 

ASD  (cols. 51-5^)  The  angular  accuracy  of  the  observation;  units  and  format 

depend  upon  the  value  of  NCODE:  for  NCODE  =  0,  1,  ^ 
or  7  ASD  is  in  seconds  of  arc  with  a  decimal  point 
implicit  between  columns  55  and  5^;  if  NCODE  =2,  5  or  8 
ASD  is  in  minutes  of  arc  with  a  decimal  point  implicit 
between  columns  52  and  55;  if  NCODE  =  5,  6  or  9  ASD  is  in 
degrees  with  a  decimal  point  implicit  between  columns 
51  and  52. 
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INEQEQ  (col. 55)  A  code  digit  indicating  the  date  of  the  equator  and 

equinox  to  which  right  ascension  and  declination  observa- 
tlona  are  referred;  0  refers  to  the  true  equator  and 
equinox  of  date,  and  1,  2,  3,  4  and  5  refer  respectively 
to  the  mean  equator  and  equinox  of  1855.0,  1875.0, 

1900.0,  I950.O  and  date.  T'^^s  digit  is  ignored  with 
observations  of  azimuth  and  elevation. 

RANGE  (cols. 56-63)  The  range  of  the  satellite,  measured  in  kilometers,  the 

decimal  point  implicit  between  columns  60  and  61;  e.g. 
587.32  km  is  punched  as  bb58732b. 

R3D  (cols. 64-68)  The  range  accuracy,  measured  in  kilometres,  a  decimal 

point  implicit  between  columns  65  and  66;  e.g.  1  .5  km  is 
punched  as  bl5bb. 

Notes ;  (i)  Other  columns  on  the  card  are  used  by  RSRS  (e.g.  "brightest 

magnitude"  in  columns  69-71 ),  but  only  those  detailed  above  are  read 
by  PROP. 


(ii)  Either  the  angular  or  the  range  part  of  the  observation  may  be 
omitted,  the  appropriate  field  being  left  blank.  The  corresponding 
"accuracy  fields"  should  also  be  left  blank,  of  course. 

(iii)  If  no  estimate  of  angular  or  range  accuracy  has  been  given  for 
an  observation,  the  appropriate  field  should  be  left  blank  when 
punching  the  card.  A  standard  value  will  be  set  from  SIGMA(4)  for 
an  angular  observation,  or  from  SIGMA(2)  for  a  range  observation 
(see  section  5.5  for  values  of  SIGMA). 

(iv)  The  zero  value  of  NCODE  (corresponding  to  a  right  ascension/ 
declination  observation,  fully  corrected  for  refraction)  was  not 
included  in  the  format  as  originally  agreed  between  the  various 
organisations  concerned.  It  was  added  at  a  late  stage  in  the 
development  of  PR0P3  and  thus  represents  an  addition  to  the  agreed 
standard  format. 


(v)  Observations  of  right  ascension  and  declination,  because  they 
are  evaluated  relative  to  a  star  background,  have  been  automatically 
corrr-cted  for  astronomical  refraction;  this  is  why  there  is  no  NCODE 
value  for  "uncorrected  right  ascension  and  declination". 
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B.7  LOOXEE  777777 

Concurrently  with  the  final  reviaion  of  the  draft  of  this  Manual  a 
seventh  observations  format  is  being  added  to  PROP.  This  provides  for  input 
of  range,  azimuth  and  elevation  for  a  named  (instead  of  numbered)  station. 
The  format  is  (2A3,  2X,  512,  f8.6,  4x,  JElJ.O)  and  the  «iuantities  on  a  card, 
in  order,  are  station  name,  date,  time,  range  (ft),  azimuth  and  elevation 
(degrees ) , 


EggtnoM  .an  1/»XEE  111111  and  UWEE  555555 

Tliese  two  formts  are  almost  Identical  and  in  Ref. 8  there  Is  In  fact  a  slnele  ronaot  for  j.A.O. 
observations.  As  stated  in  the  footnote  on  page  52  it  is  hoped  to  extend  LOOXEE  111111  to  be  in  ocoyileCs 
conformity  with  the  S.A.O.  fomsti  vhen  this  has  been  done  LOOXEE  555553  will  be  dropped.  An  initial 
extension  of  LOOXEE  111111  is  almost  complete  and  should  be  incorporated  in  PR0P3  by  the  end  of  Ilaroh,  1969* 
Ihe  format  will  be  modified  to 

FORMAT  (17,  5X,  15,  16,  212,  F6.I4,  13,  12,  F5.3,  A1 ,  212,  Ff4.2,  II,  12,  211,  7X ,  FlJi) 
and  the  following  remarks  supply  the  necessary  amending  Infomntlon  to  section  B.1. 

(1)  The  format  covers  observations  of  azlmth  and  elevation  as  well  as  right  ascension  and 
declination.  Thus  HRAH  (vSileh  now  occupies  col.3^  as  well  as  cols.35~36)  may  be  degrees  of  azlmit  1, 

Instead  of  hours  of  right  ascension,  etc. 

(2)  An  additional  variable,  (CODE,  lo  read  from  col.56.  This  may  be  0  (indicating  an  observa¬ 
tion  of  right  ascension  and  declination),  1  (for  azlnith  and  elevation,  corrected  for  refraction) 

or  3  (azlmth  and  elevation,  not  corrected). 

(3)  NTSD  (col.53)  and  NASO  (cols.5^53}  have  been  extended  to  cover  all  the  values  permitted  in 
Ref. 8,  zero  values  are  interpreted  by  0®,05  and  the  value  of  3ICMA(i*),  respectively.  NASD  is  equal  to 
the  estimated  accuracy  in  seconds  of  arc  for  values  up  to  21. 

Ihis  extension  will  enable  standard  Moonwatch  cards  to  be  read  under  LOOXEE  111111 ,  with  estlmatsd 
accuracies  used  when  available.  LOOXEE  355555  will  have  to  be  used  only  for  the  non-standard  cards  in 
which  azlmth  and  elevation  cards  are  Indicated  by  an  A  in  cols.Ul  and  5l ,  and  a  zero  in  col. 57, 

Instead  of  by  a  1  or  3  in  col. 56. 
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INTERPRETATIOW  OF  THE  SIX  MAIK  ORBITAL  EUSMEHTS 

It  may  be  important  to  the  user  of  PROP  to  know  precisely  how  to  interpret 

the  six  main  elements  of  the  PROP  model,  viz.  e  ,  i  ,  n  ,  w  ,  M  and  n 

0  0  0  0  0  0 

(where  n^  denotes  the  exclusive  element  ) .  It  has  happened  in  the  past 
that  an  organization  has  published  accurate  orbital  elements  of  a  number  of 
satellites  but  that  the  usefulness  of  these  elements  has  been  limited  by 
the  failure  of  the  organization  to  publish  definitions  of  the  elements. 

Roughly  speaking  PROP  elements  are  osculating  elements  with  short- 
periodic  pertui'bations  removed.  This  statement  may  be  made  more  precise  by 
giving  the  formulae  for  the  differences  between  PROP  elements  and  osculating 
elements,  and  this  is  done  in  the  next  paragraph.  It  must  be  stressed,  how¬ 
ever,  that  these  formulae  are  for  reference  only.  They  are  not  (apart  from 
the  formula  in  i)  used  directly  in  PROP,  which  applies  short-periodic  per¬ 
turbations  to  a  set  of  six  quantities  which  are  intermediate  between  elliptic 
elements  and  Cartesian  coordinates  -  a  set  chosen  to  make  all  the  formulae  as 
simple  as  possible  .  It  should  also  be  stressed  that  PROP  elements  have  not 
been  freed  from  long-periodic  perturbations.  In  fact  long-periodic  perturba¬ 
tions  in  PROP  (unlike  those  of  the  Smithsonian  Astrophysical  Observatory  for 
instance)  are  defined  so  as  to  be  zero  at  epoch.  This  means  that  the  only 
harmonics  involved  in  the  following  formulae  are  Jg  (for  the  short-periodic 
terms)  and  Jg  g  (for  the  along-track  correction  referred  to  in  Ref  ,2  and  in 
section  2.1  of  the  present  paper).  Also  it  is  convenient  to  regard 
as  defined  relative  to  the  true  equinox  of  date  and  hence  to  Include  a  term 

for  transformation  of  equinox  in  the  formula  connecting  n  and  fi  „  . 

o  osc 

Let 


k,  =  (jg/52)(R/p)^  and  kg  =  Jg  g  (r/p)^  /  (6  -  «)  , 

where  p  =  *1^(1  ”  ®q)  ®  modified  sidereal  angle.  To  give  first- 

order  expressions  it  is  legitimate  and  convenient  to  drop  the  zero  s\iffix 

when  writing  functions  of  e_,  i_  etc:  also  for  convenience  we  Introduce 
2  0  0 

f  =  sin  1.  Then,  if  g  is  the  longitude  of  one  extremity  of  the  major 
axis  of  the  earth’s  equator  and  v  is  the  true  anomaly  corresponding  to  the 
epoch  value  of  the  mean  anomaly,  the  required  first-oider  expressions 

are ;- 
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e  -  e  =  It.  [  6(4+e^)(2-3f )  cos  v  +  5(4+11e^)f  cos  (2«»>+v) 

OS  C  0  I 

4  (284l7e^)f  cos  (2U43v)  4  I2e(2-3f)  cos  2v 

4  6ef  {l0  cos  (3«)42v)  4  5  cos  (2i*)44v)) 

^  '  2 
4  1  cos  (2w-v)  4  cos  (3*^+5v)}  4  2e  (2-3f)  cos  3v 

4  I8ef  cos  2w  4  4e(2-3f)  {5  +  2(l-e^)/(1  4  (1-e^)^)}]  , 

i  “  i  =  4k-  sin  2i  [  3  cos  (2w42v)  4  3e  cos  (2i*>+v)  4  e  cos  (<'^45v)] 

OS  C  0  I 

^osc  ”  ^0  ”  i  [  6  (M^  -  V  -  e  sin  v)  4  3  sin  (2u»42v)  4  3e  sin  (2<»!4v) 

4  e  sin  (2u)43v)] 

4  6.1l9xlo'^  (MJD  -  35281  .9234)  4  Equation-of -Equinoxes 
[24(4-5f)(v-M^)  4  6e‘^  {  4(2-5f )  4  e^(14-l7f)}  sin  v 

Oo  C  O  I  O 

-  5e  ^  {  4f  4  (8-l5f)}  sin  (2w4v) 

4  e"^  {  (28f  -  (8-l9f)}  sin  (2w43v) 

4  l2(2-5f)  sin  2v  -  l2(2-5f)  sin  (2U42v)  4  I8f  sin  (2U44v) 

-  I8f  sin  2w  -  5ef  sin  (3*>-v)  4  2e(2-3f)  sin  3v 

5ef  sin  (2U45v)]  , 
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**oic  “  (1-e^)^  [f  {18  sin  arf  +  5e"^  (^+5*^)  sin  (2w+v) 

-  e"^  (28-e^)  sin  (2w+3v)  -  18  sin  (3J+4v) 

+  3e  sin  (2w-v)  -  Je  sin  (2w+5v)} 

-  (2-3f)  {  6e~^  (4-e^)  sin  v  +  12  sin  2v  +  2e  sin  3v}] 
-  (9  kg  n/(2*e-2fi)}  f  sin 

and 

^  -3  k-  n(l-e^)“''  [f  {l2  (2+3e^)  cos  (3»>+2v)  +  9©  (4+e^)  cos  (2w+v) 

OSC  O  I 

4  9e  (446^)  cos  (2W43v)  4  I8e^  cos  (2t»)44v) 

4  3©^  cos  (2w-v) 

TL  2  % 

4  3e'^  cos  (2w45v)  4  I8e  cos  3<>} 

f  2  2 

4  (2-3f)  { 6e  (44e  )  cos  v  4  12e  cos  2v 
4  cos  3v  +  4(243e^)}] 

It  must  be  emphasised  again  that  these  fonmilae  sire  given  only  as  a  means 

of  defining,  to  first-order,  the  PROP  elements;  the  omitted  second-order 

term  in  is  unfortunately  not  insignificant.  The  e”^  terms  in 

u  -  u  and  M  -  M  arise  because  perigee  is  not  well-defined  for  near- 
osc  osc 

circular  orbits.  This  singularity  is  unimportant  in  practice  since  PROP 
never  con^jutes  osculating  elements,  but  from  mean  elements  goes  directly  to 
Cartesian  coordinates  by  a  procedure  which  is  free  of  singularities.  The 
best  method  of  obtaining  osculating  elements  frcrni  mean  elements,  or  vice  versa, 
is  in  fact  always  to  compute  coordinates  as  an  intermediate  step^. 

Some  users  of  PROP  may  wish  to  know  how  PROP  elements  relate  to  the 
Pegasus  elements  previously  used  .  To  indicate  the  relationship  we  now 
give  a  complete  set  of  formulae  by  means  of  which  a  set  of  Pegasus  elements 
may  be  obtained  from  a  set  of  PROP  elements.  These  formulae  include 
conversions  for  'hanges  of  units. 
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Let  e,  1,  0,  h),  M,  M,  (=n),  Mo  etc.  be  a  let  of  PROP  elementii 
and  let  a',  e',  1*,  fl*,  w*,  M*,  ii|  etc.  be  the  equivalent  aet  of  Pegaius 
elements  at  the  same  epoch.  If  the  epoch  is  not  an  ascending  node  (as  it 
always  was  for  Pegasus  coiiq)uter  orbit  determinations)  the  dashed  elements  must 
be  Interpreted  not  as  true  smoothed  elements,  as  originally  defined  by 
Merson,  but  as  mean  smoothed  elements  (related  to  tme  smoothed  elements  in 
the  same  way  that  PROP  elements  are  to  osculating  elements). 

Taking  angles  in  degrees  we  require  a  constant  k  (=  57.2957795)  for 
conversion  from  radians.  We  also  require  conversions  in  time  since  is  in 

"/day/day  but  n^  is  in  •/(!  00  days)  .  Let  k^  be  the  same  constant  as 
before.  As  in  the  previous  formulae  zero  suffices  are  omitted  when  possible. 

Then  a  complete  set  of  first-order  formulae  is:- 
a*  =  (86400  k/M|  )^}^  ,  where  =  398602  , 

^  =  (J2/^2)(R  /  a'(l-e^))^  , 

f  =  sin^  i  , 

e*  ^  (2-3f)(1-e^)  e"^  {(1+e^)  -  (1-e^)^  , 

i'  -  i^  -  4kk|  sin  2i  (3  +  4  e  cos  w)  , 

-  -2k  tan  ^  {  (1-e)^  (1+e)”^  tfji  ,  with  (E^  +  w|  <  I80  , 

M  =  E  -  k  e  sin  E  > 

00  0  0 

fi*  +  l6k|  cos  i  { 3(<*’  +  M^^  )  +  4k  e  sin  w} 

+  5.506  X  lo"^  (MJD  -  35281  .9254)  +  e/240  , 

where  E  is  the  equation  of  equinoxes 

w»  ^  4  4k^  {k  e^  (2-3f)  sin  2w  -  2(4-5f)(3u  4  3  4  2k  e  sin  w)}  , 

*If  e  is  small  it  is  better  to  replace  e”  {  (l4e)^  -  (l-e^)^  by  the 

2  4 

non-singular  truncated  expression  ^  (5  +ie  ). 


07 
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^5 
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52 
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53 
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54 
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45 
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57 
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10, 

62 
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51, 

57 
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64 
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6, 

10, 

16, 

17, 

2^ 

56 
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29, 

56 
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15, 

?L 

22, 

56, 

44 

equator 

9, 

26, 

55. 

60 
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9, 

26, 

55, 

60 

failure  of  PROP  (standa.d) 

16, 

17, 

45 
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16, 

17, 

45 
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light-time  correction 
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low-perigee  satellites 
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mode  of  refinement 
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Moonwatch  observations 
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observing-stations 
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of  observations 

observation-cards 
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output,  printed 

punched-card 

parameters  (orbital) 
partial  derivatives 
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prediction  of  elements 
PREP 

presentation  of  results 
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Fig.  I  Calling  structure  tor  PROPS  program  units 


Fig. 2 
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Fig. 2  Flow  chart  for  Prop  3  main  program 


Fig.  3  Voriolion  of  density  scale  height,  H,with 
oltitude  (S  =  10-7  cm  solar  radiation  in  Wm  ^  Hz"') 
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Fig. 5  Output  for  typical  PROP  run 
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